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ABSTRACT — — 

DESAP 2 is a finite element program for computer-automated, minimum weight design of 
elastic structures with constraints on stresses (including local instability criteria) and buckling 
loads. No limits are placed on the number of load conditions for stress-constrained design, but 
only one of these load conditions can be chosen as the potential buckling load. A substantial por- 
tion of DESAP 2., particularly the analysis of the prebuckling state, is derived from the SOLID 
SAP finite element program developed at the University of California, Berkeley. The stress- 
constrained design is based on the classical stress ratio method, which drives the design toward 
a fully stressed state. The constraints on the buckling load are handled by solving the appropri- 
ate optimality criterion by successive iterations. During each iteration, the element sizes deter- 
mined by the stress ratio method are used as the minimum size constraints. The element 
subroutines have been organized in a manner that permits the user to make additions and changes 
with a minimal programming effort. Consequently, DESAP 2 can readily be changed into a 
special-purpose program to handle the user f s specific design requirements and failure criteria. 

DESAP 2 is a companion program of DESAP i: "A Structural Design Program with Stress 
and Displacement Constraints. M With the exception of a few cards the same input data deck can 
be used for both programs. 
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K. PREAMBLE 


This collection of sample problems is a supplement to Volume 1 
of DESAP 2: ’Theoretical and User's Manual". In making up this 

volume, our aim was to find examples that would best serve the 
following functions: 

1) Illustrate and supplement the input instructions of Volume 1, 
and to familiarize the user with the output. 

2) Explain, with examples, special problem areas and peculiar- 
ities that may arise in the use of the program. 

3) Provide example problems that may be used for debugging the 
program during installation on a new computer system. 

4) Compare the results of DESAP 2 against solutions obtained by 
other means, whenever possible. 

Although DESAP 2 is designed primarily for the use of large 
structures, the stated purpose of the sample problems is clearly best 
fulfilled by small, simple examples that do not necessarily represent 
realistic design situations. Consequently, the problems appearing 
in this volume should be viewed strictly as tools of instruction, 
which in no way reflect the ultimate capabilities of the program. 

Because our experience with the program is rather limited at 
this time, the example problems may well have overlooked some trouble- 
some aspects of design, or even deficiencies in the program itself. 

The extensive computer output from each design cycle is, however, a 
powerful diagnostic tool that should enable the user to pinpoint the 
difficulty and make the appropriate correction. 


K. 2 


An example problem is given for each element type presently used 
in the program. Each problem contains a complete description of the 
input data, including an echo of the input cards, and the computer 
printout of the input information. In order to reduce the bulk of 
the report, only a partial listing of the computer output is duplicated, 
containing the initial and the final designs. The complete history of 
a design is usually summarized by tabulating the design variables. 

In compiling the sample problems, we were seriously handicapped 
by a lack of adequately documented optimal design problems in existing 
literature. For this reason, a one-to-one comparison of the results 
of DESAP 2 with independently obtained solutions is lacking in some of 
the problems. 

As a final note, we would like to remind the user again that 
DESAP 2 is oriented towards large problems. Mainly due to an extensive 
use of auxiliary storage devices and other core-saving features, the 
program is not efficient for small structures as used for the sample 
problems. Consequently, the computer times for these problems are 
not expected to be competitive with runs obtairied from programs 
especially designed for structures of limited size. 
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L, BAR ELEMENTS 


L.l Three-Panel Truss 



Layout of Truss Showing Element and Node Numbers. 


The truss in Fig. L.1.1 is to be optimized for the single load 
condition shown. The cross-sectional area of each bar is taken as 
an independent design variable. The data used in the design is: 

E = 10 psi (Young’s modulus), 

a* = a* = 25 psi (allowable normal stress), 

A* = 0.1 sq. in. for all elements (min. allowable 
cross-sectional areas) , 

A = 1.0 sq. in. for all elements (initial cross- 
sectional areas). 


L, 1 . 2 


p* = 1.0 (lower bound on the critical load parameter, i.e. 
the factor of safety against buckling). 

Local buckling of the elements is not to be included as a design 
criterion. Because the two construction codes of the bar element 
differ only in redesign with respect to local buckling, either of the 
codes may be used (we chose Construction Code No. 1). 

The complete design history of the problem is given on the 
Computer Printout sheets. The optimal design was reached in two re- 
designs; it is governed entirely by the buckling constraint. Because 
stress constraints are inactive and the prebuckling state is statically 
determinate, the final design is a true global optimal design. 

The buckling load of the initial design (see Analysis of Design 
No. 0) of 2.2654 agress exactly with the analytical solution of Ref. [13], 
p. 148. Unfortunately, the cross-sectional areas of optimal design 
cannot be checked due to lack of an independently obtained solution. 
Special notes on input -output : 

1) Since the stiffness matrix of each element of the structure has the 
form [IC] = [k^jA?, n = 1, uniform scaling is an exact operation. 
Consequently, KSCALE = 1 (=n) was specified in the Design Control 
Data. 

2) Automatic generation of boundary condition codes was employed to 
suppress the z-displacement and the rotations at all the nodes 
(see Nodal Point Input Data and Generated Nodal Data). 

3) Local buckling of elements was eliminated as a design consideration 

by leaving the moments of inertia blank on the geometric property 

6 4 

card. The blanks are replaced by the computer with 10 in 
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(see Geometric Property Cards), so that the Euler buckling load 
of each element is too high to play a role in the redesign. 

4) MODEIN = 1 in the Buckling Control Data means that the initial 
mode shape (coordinate vector) is read Xn with data cards, 

5) NMODE = 1 in the Buckling Control Data specifies that only one 
buckling mode is to be considered in the design. 

6) The parameter INDET = 1 in the Buckling Control Data signifies 
that the structure is 'statically determinate. Therefore, the pre- 
buckling stresses and the geometric stiffness matrix of the struc- 
ture have to be calculated only once in the whole redesign process. 

7) A "normal" value of the relaxation parameter a = n/(n+l) = 0.5 
(note that n = 1) was employed in the buckling-constrained redesign 
(see ALPA in Buckling Control Data) . 

8) The program was terminated when the Optimality Index of each ele- 
ment became sufficiently close to one (see Evaluation of Design 
No. 2). 


18400 1 2 34 5ft 7 8') A ] 2 34 56789 Ml 2 3456 7S9CI 2 34 5* 7H9I) 1 2345 67R9KI 234 56 789 FI 2 3456 78901 2 34 5 6 78911 
13450 TURK EE PAN El. TRUSS BUCKLING ANALYSIS 
18500 5117 

1H550 10 l 0.025 I l 

1 S 600 1 l 1 “1 - l - l -1 


18630 

2 

20. 


187UO 

3 

10.0 

10 

1 8 7 5 0 

4 1 

40 . 


1 8800 

5 1 

1 1 1 30. 

10 


18850 171111 

1 H 9 0 0 1 l 0-1 

18930 10. 25. 

1 9 UOU 1 1. 

19 0 5 0 
19100 
19150 
192 00 

19250 1 1 2 l 1 1 

19300 2 2 A 1 12 

19350 3 3 5 l l 3 

19400 413114 

19450 523115 

19500 625116 

19550 745117 

1 9600 

19650 1.0 1 l l I 0.5 0.8 

19700 3 1 1. 

19750 4 1 -1. 

19800 5 l -I, 

19850 

19900 0.5 1.0 -1.0 1.0 -1.0 1.0 1.0 

19950 3 1.0 O.l 

20000 7 2.0 O.l 

200 50 1 2 34 5678 9 A123 4 56789R12 3456789 C 1274 5 6789D I 234567 S9C12 34 5G789FI 2345 6 7 89C I 2345678 911 

2 010 0 

Echo of Input Cards 




THFEEE PANEL TRUSS BUCKLING ANALYSIS 


NUMBER OP NODAL POINTS - 5 
NUMBER OP ELEMENT TTPF.S = 1 
NUMBER OF LOAD CASES = 1 
NUMBER OF DES. VARIABLES = 7 


DESIGN CONTROL DATA 

NCTCL - 10 

KSCA IE= 1 

DELTA = 0.2500E-01 

EPSIL = 0 . 1000 E 00 

LBUCK - 1 


NODAL POINT INPUT DATA 


NODE 

BOUNDARY 

CONDITION 

CODES 

/ 

NODAL POINT 

COORDINATES 

NUMBER 

X 

Y 

Z 

XX 

YY 

zz 

X 

T 

1 

1 

1 


-1 

-1 

-1 

0.000 

0.000 

2 

0 

0 

0 

0 

0 

0 

20. 0C0 

0.000 

3 

0 

0 

0 

0 

0 

0 

10.000 

10.000 

U 

0 

1 

0 

0 

0 

0 

40.000 

0.000 

5 

0 

0 

1 

1 

1 

1 

30,000 

10.000 


GENERATED 

NODAL 

DATA 






NODE 

BOUNDARY 

CONDITION 

CODES 

/ 

NOC AL POINT 

COORDINATES 

NUMBER 

X 

Y 

Z 

XX 

YY 

ZZ 

X 

Y 

1 

1 

1 

-1 

-1 

-1 

-1 

C.000 

0.000 

2 

0 

0 

-1 

-1 

-1 


2C.OOO 

0.000 

3 

0 

0 

-1 

-1 

-1 


1C. 000 

10.000 

4 

0 

1 


-1 

-1 

-1 

40.000 

0.000 

5 

0 

0 

1 

1 

1 

1 

30.000 

10.000 


EQUATION 

NUMBERS 





N 

X 

T 

Z 

XX 

YY 

zz 

1 

0 

0 

0 

0 

0 

0 

2 

1 

2 

0 

0 

0 

0 

3 

3 

4 

0 

0 

0 

0 

4 

5 

0 

0 

0 

0 

0 

5 

6 

7 

0 

0 

0 

0 


Computer Printout 


/ 

Z T 


0.000 

0 

0.000 

0.000 

0 

0.000 

0.000 

0 

0.000 

0.000 

0 

0.000 

0.000 

0 

0.000 


/ 

Z T 


0.000 

0 , 

,000 

0.000 

0 , 

,000 

0.000 

0 , 

,000 

0.000 

0 . 

,000 

0.000 

0 , 

,000 


In 



NUMBER OF TRUSS EL EHENTS = 7 
CONSTR UCTIOH CODE = 1 
NUMBER OF MATERIALS =* 1 
HUMBER OF TEMPS FOR WHICH MATL PROPS GIVEN= 1 
NUMBER OF DIFFERENT GEOMETRIES PROPS GIVEN= 1 


MATERIAL PROPERTY CARDS 


MATERIAL NUMBER SPECIFIC YOUNGS COEFFT OF /—ALLOWABLE STRESSES — / 

NUMBER OF TEMPS WEIGHT TEMP MODULUS THERM EXPAH TENSION COMPRESSION 

1 1 0.1000E 00 0.0000E 00 0.1000E 02 0.0000R 00 0.2500E 02 0.2500E 02 


GEOMETRIC PROPERTY CARDS 

GEOMETRY X-SECT /—MOMENTS OP INERTIA—/ 

NUMBER AREA YY ZZ 

1 0.1000D 01 0.1000E 07 0. 1000E 07 


ELEMENT LOAD MULTIPLIERS 




A 


B 


C 


D 

X-DIR 

0.000000D 

00 

0. OOOOOOD 

00 

0. OOOOOOD 

00 

0. OOOOOOD 

00 

Y-DIF 

0.000000D 

00 

0. OOOOOOD 

00 

0. OOOOOOD 

00 

0. OOOOOOD 

00 

Z-DIR 

0.000000D 

00 

0. OOOOOOD 

00 

0. OOOOOOD 

00 

0. OOOOOOD 

00 

TEMP 

0. OOOOOOD 

00 

0. OOOOOOD 

00 

0. OOOOOOD 

00 

0. OOOOOOD 

00 


PROCESSED 

ELEMENT 

DATA 













ELEMENT / 

-NODE NOS-/ 

/--ELEMENT 

ID 

NOS-/ 

DESIGN VAR 

RFFEPENCE 

END FIXTTY 

COEFFICIENTS 

BAND 

NUMBER 

I 

J 

MATL GEOMY 

D VAR 

FRACTION 

TEMP 


YY 


7.Z 


WIDTH 

1 

1 

2 

1 

1 

1 

0. 1000E 

01 

0.0000D 

00 

0.1000D 

01 

o.iooon 

01 

2 

2 

2 

4 

1 

1 

2 

O.IOCOE 

01 

O.OOOOD 

CO 

0.1000D 

01 

O.IOCOD 

01 

S 

3 

3 

5 

1 

1 

3 

0. 1 0 COE 

01 

O.OOOOD 

00 

0.10C0D 

01 

0.1000D 

01 

5 

4 

1 

3 

1 

1 

4 

0. 1000E 

01 

O.OOOOD 

00 

0.10C0D 

01 

0.1000D 

01 

2 

5 

2 

3 

1 

1 

5 

0. 1000E 

01 

O.OOOOD 

00 

0.1000D 

01 

O.IOCOD 

01 

4 

6 

2 

5 

1 

1 

6 

0. 1000 E 

01 

O.OOOOD 

00 

0. 1000D 

01 

0.1000D 

01 

7 

7 

4 

5 

1 

1 

7 

C.1000E 

01 

O.OOOOD 

00 

O.IOCOD 

01 

0.1000D 

01 

3 


STRUCTURE 
LOAD CASE 

1 


STRUCTURE LOAD MULTIPLIERS 
A B C D 

0.000 0.000 O.OCO 0.000 


BUCKLING CONTROL DATA 


COEFFT 

HODEIN 

NMODE 

INDET 

NVEC 

ALPA 

OMEGA 


1.00000 

1 

1 

1 

1 

0.50000 

0.80000 


tr 


ON 



NODAL POINT LOADS 


o.ooon oo 

O.OOOD 00 

o.ooon oo 


DESIGN VARIABLE INPOT DATA 
DESIGN 

VARIABLE INITIAL BIN ALLOWABLE 

NUMBEP VALUE VALOE 

1 0. 1000E 01 0. 1000E 00 

2 0.1000E 01 0.1000E 00 

3 0.1000E 01 0. 1000E 00 

4 0.2000E 01 0.1000E 00 

5 0.2000E 01 0. 1000B 00 

6 0.2000E 01 0.1000E 00 

7 0.2000E 01 0.100QE 00 


TOTAL HUMBER OP EQUATIONS = 7 
BANDWIDTH = 7 
NUMBER OF EQUATIONS IN A BLOCK = 7 
NUMBER OF BLOCKS * 1 


NODE 

LOAD 



APPLIED 

LOADS 


NO. 

CASE 

FX 

RY 

RZ 


3 

1 

0. 100 D 

01 

O.OOOD 00 

O.OOOD 

00 

4 

1 

-0. 100D 

01 

O.OOOD 00 

O.OOOD 

00 

5 

1 

-0. 100 D 

01 

O.OOOD 00 

O.OOOD 

00 


MY MZ 

c.oood oo o.oood no 
O.OOOD 00 O.OOOD 00 
0. ooon CO 0 .OOOd 00 


t- 


■^1 



****************************** 
AS ALTS IS OF DESIGN Nil MBER 0 

****************************** 


HODAL DISPLACEMENT S AND ROTATIONS 


NODE 

NO, 

load 

CASE 

X 

Y 

Z 

XX 

YT 

ZZ 

5 

1 

-3.000E 00 

1.000E 00 

0 * OOOE- 0 1 

O.OOOOE-OI 

C.0000F-01 

O.OOOOE-OI 

4 

1 

-4. OOOE 00 

O.OOOE-OI 

0.0C0E-01 

0 • OOOOE-O 1 

0.0000E-01 

O.OOOOE-OI 

3 

1 

-1.000E 00 

1.000E 00 

O.OOOE-OI 

O.OOnoE-OI 

O.OOOOE-OI 

O.OOOOE-OI 

2 

1 

-2. OOOE 00 

4 . 330E-15 

O.OOOE-OI 

O.OOOOE-OI 

0 . OOOOE-O 1 

O.OOOOE-OI 

1 

1 

0 • 000E-01 

0. 000E-01 

0.000E-01 

0.0000E-C1 

O.OOOOE-OI 

0 . 00 00E-01 

VALUES OF 

DESIGN VARIABLES 







1 

2 

3 

4 

5 

6 


0 0.1000E 01 0. 1000 E 01 0. 1000E 01 C.2000F 01 0.2000E 01 0.2000E 01 0.2000E 01 


ANALYSIS OF TRUSS ELEMENTS , CONSTPN CODE= 1 
ELEMENT X-SECT AREA LOAD COND AXIAL FORCE 


1 0.1000E 01 

2 0.1000E 01 

3 0.1000E 01 

4 0.2000E 01 

5 0.2000E 01 

6 0.2000E 01 

7 0.2000E 01 


-0.1000E 01 
-0.1000E 01 
-0.1000E 01 
O.5969E-07 
0 . 9537E- 06 
-0.9537E-C6 
0.9537E-06 


BUCKLING LOAD PARAMETERS 
0.22654D 01 


BUCKLING MODE SHAPES 


NODE MODE 
NO. SHAPE 

X 

T 

7 , 

XX 

YY 

Z7, 

5 1 

5. 895E-0B 

-1.939E 00 

O.OOOE-OI 

0 . 000 n F- Cl 

o.oooon-oi 

O.OOOOE-OI 

4 1 

1.433E 00 

0 . OOOE- 0 1 

0 . OOOE- 0 1 

0 . OOOOE-O 1 

0.0000E-01 

0 . 00 00 F- 0 1 

3 1 

1.433E 00 

-1.939E 00 

O.OOOF-OI 

0 .0000 E-01 

O.OOOOF-C1 

0 . CO 00 E- 0 1 

2 1 

7. 16UE-01 

-3. 162E 00 

0.Q00E-01 

0 ,OOOOE-O1 

c.oproF-ci 

0 . OOOOE-O 1 



1 


1 O.OOOF-OI 0. 000E-Q1 0.000E-01 0.0000E-01 0.0000E-C1 O.OOOOE-OI 


***** ***************** ********* 
EVALUATION OF DESIGN NUMBER 0 
******************************* 


STRESS RATIO LOAD COND DES VARIABLE 
MAX 0. 1000E 00 0 1 

MIN 0. 5000E-0 1 0 4 

BAX BOCK RATIOS LOAD COND 

0.4414F 00 1 


UNIFORM SCALING OPERATION FOLLOWS 


SCALE FACTOR IS 0.441AKD DETERMINED BY BUCKLING CONSTRAINTS 


DESIGN VARIABLES OF SCALED (CRITICAL) DESIGN ARE 


VALUES OF DESIGN VARIABLES 

1234567*0 10 

0 Q.4414E 00 0. 44 14 E 00 0.4414E 00 9.882RE 00 0.P828E 00 0.8R2RE 00 0.8828E 00 


STRUCTURAL WEIGHT^ 0.7641E 01 


REDESIGN OPERATION FOLLOWS 


OPTIMALITY INDEX OF DESIGN VARIABLES FOR BUCKLING CONSTRAINTS 
DV NO ACT/PAS INDEX 

1 ACT 0.45474E 00 

2 ACT 0. 45474E 00 

3 ACT 0. 18190 E 01 

4 ACT 0.22737B 00 

5 ACT 0.22737E 00 

6 ACT 0.2273 7E 00 

7 ACT 0. 227375 00 

NO. OF ACTIVE BUCKLING CONSTRAINTS ARE 1 




to 


i 


****************************** 

ANALYSIS OF DESIGN NUMBER 1 

****************************** 


NODAL DISPLACEMENTS AND ROTATIONS 


NODE 

NO. 

LOAD 

CASE 

X 

Y 

7 

XX 

YY 

7.Z 

5 

1 

-7.R36E 00 

4.622E 00 

O.OOOE-OI 

O.OOOOE-OI 

0.0000F-C1 

O.OOOOE-OI 

4 

1 

-1.246E 01 

0.000B-01 

C • 000E-C1 

0 • OOOOE-O 1 

O.OOOOE-OI 

0.0000E-C1 

3 

1 

-4.622E 00 

4.622E 00 

O.OOOE-OI 

O.OOOOE-OI 

O.OOOOE-OI 

0 . OOCQE-C 1 

2 

1 

-6.229E 00 

3.015B 00 

O.OOOE-OI 

0 . OOOOE-O 1 

0. OOOOE-O 1 

0. OOOOE-C1 

1 

1 

0.000E-01 

O.OOOE-OI 

O.OOOE-OI 

O.OOOOE-OI 

0.0000F-01 

0.00C0E-01 

VALUES OF 

DESIGN VARIABLES 







1 

2 

3 

4 

5 

6 


0 0.3211E 00 0.3211M 00 0.6222E 00 0.541BE 00 0.5418E 00 0.541BE 00 0.5418E 00 


ANALYSIS OF TRUSS ELEMENTS , CONSTFN CODE= 1 


ELEMENT 7-SECT AREA LOAD COND AXIAL FORCE 


1 0. 321 1 E 00 

2 0.3211E 00 

3 0.6222E 00 

4 0.5419E 00 

5 0 . 54 1 BE 00 

6 0.5418E 00 

7 0.5418E 00 


-0.1000E 01 
-C.1000E 01 
-0.1000E 01 
0.2347E-06 
0.1729E-C5 
-C. 2603E-C6 
0. 9689E-06 


BUCKLING LOAD PARAMETERS 
0.86583D 00 


BUCKLING MODE SHAPES 


NODE MODE 
NO. SHAPE 

X 


Y 

7 

XX 

YY 

77 

5 

1 

-4.127E-01 

2.00BE 

00 

3.000E-01 

O.onoOE-OI 

C.000CF-01 

O.OOOOE-OI 

4 

1 

-1.706E 00 

0.000E- 

•01 

O.OOOE-OI 

O.OOOOF-OI 

O.OOOOE-OI 

0 . 0000 E- 0 1 

3 

1 

-1.293E 00 

2.008E 

00 

O.OOOE-OI 

0.00C0E-01 

0.0090F-01 

O.OCOOE-C1 

2 

1 

-8.52BE-01 

3. 162E 

00 

C.900E-01 

9.0000F-01 

0. OOOOE-O 1 

0 . OOOOE-O 1 


10 


L.1.10 


1 


1 0 . OOOE-O 1 0.000E-01 


C.000F-01 O.OOOOE-OI 


O.OOOOE-OI 0 . COOOE-OI 


******************************* 
EVALUATION OF DESIGN NUWBER 1 
******************************* 


WAX 

WIN 


STPESS RATIO 
0.3115E 00 
0.1607E 00 


LOAD COND PES VARIABLE 

0 2 

0 3 


WAX BUCK RATIOS LOAD COND 


0.1155E 01 1 


TINIFOFB SCALING OPERATION FOLLOWS 


SCALE FACTOR IS 1.155AND DETERMINED BT BUCKLING CONSTRAINTS 


DESIGN 

VARIABLES 

OF SCALED 

(CRITICAL) DESIGN ARE 



VALITES 

OF DESIGN 

VARIABLES 





1 

2 

3 4 

5 6 

7 

0 

0.3708E 00 

( 0. 3703 E 

00 0.7186E 00 0.6257E 00 

0.6257E 00 0.6257E 00 

0.6257E 00 


STRUCTURAL WEIGHT^ 0.6460E 01 
REDESIGN OPERATION FOLLOWS 


OPTIMALITY INDEX OF DESIGN VARIABLES FOP BUCKLING CONSTRAINTS 
DV NO ACT/PAS INDEX 


1 

ACT 

0.11317E 01 

2 

ACT 

0.11J17B 01 

3 

ACT 

O.12055E 01 

4 

ACT 

0.79490E 00 

5 

ACT 

0.79490E 00 

6 

ACT 

0.79490E 00 

7 

ACT 

0. 79490E 00 


ABE 


NO. OF ACTIVE BUCKLING CONSTRAINTS 


1 



****************************** 
ANALYSIS OP DESIGN NOHBFR 2 
****************************** 


NODAL DISPLACEMENTS AND ROTATIONS 


NODE LOAD X 

NO. CASE 

Y 

55 

XX 

IT 

zz 

5 1 -6.322E 00 

3.79BE 00 

O.OOOE-OI 

O.OOOOE-OI 

O.OOOOF-OI 

O.OOOCE-OI 

U 1 -1.012E 01 

0. OOOE-O 1 

3 • OOOE-0 1 

O.OOOOE-OI 

O.OOOOE-OI 

O.OOOOE-OI 

3 1 -3.798E 00 

3.798E 00 

0 .OOOE-O 1 

0 . OCOOE-0 1 

O.OOOOF-OI 

0 . 0000 E- 01 

2 1 -5.060E 00 

2.538E 00 

0.000B-01 

O.OOOOE-OI 

O.OOOOE-OI 

O.OOOOE-OI 

1 1 0. OOOE-O 1 

O.OOOE-OI 

O.OOOE-OI 

0 . OOOOF-O 1 

0. OOOOE-0 1 

O.OOOOE-OI 

VALDES OF DESIGN VARIABLES 





1 

2 

3 

4 

5 

6 

0 0.3952E 00 0.3952B 00 0.7924E 00 0.5616E 00 0.' 

5616E 00 0.5616E 00 0. 

ANALYSIS OF TRUSS ELEMENTS, CONSTPN CODE= 1 




ELEMENT X-SECT AREA 

LOAD COND 

AXIAL FORCE 




1 0.3952E 00 

2 0.3952E 00 

3 0 .79 24 E 00 

4 0.5616E 00 

5 0.5616B 00 

6 0.5616E 00 

7 0.5616E 00 

1 

1 

1 

1 

1 

1 

1 

-0.1000E 01 
-0.1000E 01 
-0.1000E 01 
0.74 12E-06 
0. 2309E-05 
0. 8326E-07 
0.5191E-06 




BUCKLING LOAD PARAMETERS 





0. 99 10 2D 00 






BUCKLING MODE SHAPES 






NODE MODE X 

NO. SHAPE 

Y 

7j 

XX 

YY 

Z7. 

5 1 3.974E-01 

-1.978E 00 

O.OOOE-OI 

O.COOCE-OI 

C . 0000 F- 01 

O.OOOOE-CI 

4 1 1.586E 00 

0 • OOOE-O 1 

0 .OOOE-O 1 

0.000CE-01 

O.OOOOE-CI 

O.OOOOE-CI 

3 1 1.188E 00 

-1.978H 00 

0 .POOF -01 

O.OOOQF-OI 

0.0000F-C1 

0 . CO 00 E- 0 1 

2 1 7.929E-01 

-3. 162E 00 

O.OOOE-OI 

O.OOOOF-OI 

O.OOCOE-C1 

O.COOCF-OI 


ZX‘l'1 



1 


1 0. OOOE-O 1 O.OOOB-OI O.OOOE-OI 0.0000E-01 


C.OOOOE-OI 


O.OOCOE-OI 


EVALUATION OP DESIGN NUMBER 2 

***************************+*** 


RAX 

niN 

STRESS PATIO 
0.2530E 00 
0.1262E 00 

LOAD CCHD 
0 
0 

DES V ARI ABLE 
2 
3 


RAX BUCK RATIOS 

LOAD CCND 



0.1009E 01 

1 


DESIGN 

IS CRITICAL 




STRUCTURAL WEIGHT* 0.6343E 01 


REDESIGN OPERATION FOLLOWS 


OPTIHALITT INDEX OF DESIGN VARIABLES FOR BUCKLING CONSTRAINTS' 


DV NO ACT/PAS 

1 ACT 

2 ACT 

3 ACT 

*1 ACT 

5 ACT 

6 ACT 

7 ACT 


INDEX 

0. 1024 IE 0? 
0. 1024 1 E 01 
0. 10191E 01 
0. 10146E 01 
0.10146E 01 
0. 10146E 01 
0.10146E 01 


NO. OF ACTIVE BUCKLING CONSTRAINTS ARE 1 
BUCKLING - CRITICAL DESIGN HAS CONVERGED 


L.1.13 
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L. 2 Ten-Panel Truss 



The truss in Fig. L.2.1, which has been treated in existing 
literature [14] , gives an opportunity to check the results of 
DESAP 2 against an independently obtained solution. Both designs 
use the following data: 

E = 10 x 10^ psi [Young ! s modulus), 
a* = a* = 20,000 psi [allowable stress), 

A = 2.0316 sq. in. for all members (initial cross-sectional 
area) , 

A* = 0.0544 sq. in. for all members (minimum allowable 
cross-sectional area). 

Local buckling of the elements was not taken into consideration. 


L.2.2 


The designs obtained from DESAP 2 and Ref. [14] after five 
redesign cycles are listed in Table L.2.1; the correlation is excellent. 
It should be noted that the design is symmetric, as it should be, i.e. 
the two horizontal (chord) members in each panel are equal, and so 
are the two diagonal members. 

This problem differs significantly from the three-panel truss in 
Sec. L. 1 . Firstly, the prebuckling state is not statically determinate 
in the current problem, which means that a prebuckling analysis must 
be carried out, and the geometric stiffness matrix of the structure 
recomputed prior to each design cycle. Secondly, the final design is 
governed by stress and buckling constraints simultaneously, whereas 
only the buckling constraint was active for the three-panel truss. 



L.2.3 


Panel 

number 

Cross-sectional areas fsq. in. 

) 

DESAP 2 

Reference [14] 

Hori z . 
membs . 

Diag. 
membs . 

Hori z . 
membs . 

Diag. 
membs . 

1 

2.713 

0.054 

2.721 

0.054 

2 

2.654 

0.054 

2.655 

0.054 

3 

2.536 

0.075 

2.532 

0.076 

4 

2.362 

0.103 

2.357 

0.102 

5 

2.134 

0.130 

2.132 

0.128 

6 

1.855 

0.155 

1.857 

0.153 

7 

1.527 

0.179 

1.531 

0.179 

8 

1.152 

0.203 

1.155 

0.204 


0.728 

0.229 

0.729 

0.230 

mm 

0.251 

0.256 

0.251 

0.256 

Weight 

466.4 

■ lb. 

466.6 

lb. 


The cross-sectional areas of all vertical members are at 
the minimum allowable value of 0.0544 sq. in. 


Table L.2.1 

Comparison Between the Results of DESAP 2 and Reference [14] 
After Five Design Cycles. 
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M. BEAM ELEMENTS 


M. 1 Three-Member Frame with Elastic Supports 


D 

O 



Layout of Frame Showing Element and Node Numbers 


The spring -supported frame in Fig. M.1.1, which has been 
extensively discussed in Ref. [4], is an example of a design problem 
that is governed by two buckling modes simultaneously. 

The structural members consist of thin-walled beams. Only the 
wall thickness is to be varied during the design, which means that 
Construction Code No. 1 must be used. Equal size constraints are 



M.1.2 


used to reduce the number of design variables to two: the cross- 

sectional area of the columns (A^), and the cross-sectional area 
of the beam (A 2 ). The elastic supports are modelled by boundary 
elements . 

The following data is used in the design: 

E = 30 x 10^ psi (Young 1 s modulus), 

a* = a* = 50,000 psi (allowable stress), 
x c 

p =1,0 lb./cu. in. (specific weight), 

p* = 1.0 (min. allowable critical load parameter). 

The geometric properties of the elements are summarized below. 


Design 

Variable 


Initial 

Values 

Elements 

A(in. 

I z (in. 4 ) 

i 

1-8 

1.50 

1.50 

2 

9 

0.15 

0.15 


No minimum size constraints were placed on the design variables. 

The buckling of this simple frame can be treated analytically 
for arbitrary values of the two design variables A^ and The 

results of the analysis have been plotted in the form of "design 
space" in Fig. M.1.2. Each point in the design space has as its 
coordinates A^ and and represents, therefore, a specific design. 
The constraints divide the design space into two regions: the 

feasible region consists of designs that do not violate any con- 
straints, whereas the infeasible region violates at least one 
constraint. The boundary between the two regions represents critical 
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designs. The optimal design is that point in the feasible region that 

has the lowest weight. By inspection of Fig. M.1.2 it can be seen that 

this point lies at the intersection of the two buckling constraint 

lines P = 10k and P = 10k, where P and P are the buckling loads 
s a s a 

associated with the symmetric and asymmetric inodes, respectively. 

The designs obtained by DESAP 2 are also shown in Fig. M.1.2. As 
can be seen, the optimal design was reached in four redesigns, including 
two uniform scaling operations. The history of the design process is ' 
also given in numerical form in Table M.1.4. The stress constraints 
were not active at any stage of the design. 



Design Number 

0 

1 

2 

3 

4 

Aj (in 2 ) 

1.500 

1.318 

0.858 

0.925 

0.841 

A 2 (in 2 ) 

0.150 

0.132 

0.231 

0.249 

0.289 

P s (kips) 

16-44 

14.44 

9.97 

10.74 

10.00 

P a (kips) 

10.47 

9.94 

9.65 

9.98 

9.97 

Wt. (lb) 

N.C. * 

656.2 

N.C. 

488.7 

455.8 


^Denotes that design is not critical. 

Table M.1.4 

Design History of Frame 








M.1.5 


Special notes on input -output : 

1) Although the stiffness matrix of each element has the form 

[K^] = the presence of the elastic supports means that 

uniform scaling is not an exact operation. Consequently, 

KSCALE = 0 is used in Design Control Data, which causes each 
scaled design to be analyzed and evaluated. 

2) As each element bends only about the local z-axis, the properties 
about the two other axes are left blank on the Geometric Property 
Cards. 

3) All section moduli have been left blank on the Geometric Property 
Cards. Consequently, the bending stresses are not calculated in 
the prebuckling state (they would be very small anyway). 

4) The use of node Nos. 1 and 12 as the "third 11 node on the element 
cards (see Processed Element Data) specifies local y-axis to lie 
in the global x-y plane. 

5) NMODE = 2 in Buckling Control Data requests that the redesign be 
carried out with respect to two modes simultaneously. 

6) The prebuckling state is statically determinate (bending is 
negligible), which is specified by INDET = 1 in the Buckling Control 
Data. Consequently, the geometric stiffness matrix of the structure 
is assembled only once in the entire design process. 
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T H PE EE WEHBEP FRAME WITH ELASTIC SUPPORTS 


RUCKLING DESIGN 


NUMBER OF NODAL POINTS = 12 
NUMBER OF ELEMENT TYPES * 2 
NUMBER OF LOAD CASFS = 1 
NnMBFF OF DES. VARIABLES = 2 


DESIGN CONTROL DATA 


NCYCL = 
NSC A I F“ 
DELTA = 
EPSTL = 
LDDCK = 


10 

0 

0.2500E-01 
0.1000E 00 
1 


NODAL POINT INPUT DATA 


NODE 

BOUNDARY 

CONDITION 

CODES 

/-- 

NODAL POINT 

COORDINATES-- 


•/ 



NIJNDSP 

X 

Y 

Z 

XX 

YT 

77, 

X 

Y 


Z 


T 

1 

1 

1 


-1 

-1 

1 

0.0 

0.0 

0.0 


0 

0.0 

2 

1 

1 

0 

0 

0 

1 

180.000 

0.0 

0.0 


0 

0.0 

3 

0 

0 

0 

c 

0 

0 

0.0 

60.000 

0.0 


0 

0.0 

9 

e 

0 

0 

0 

n 

0 

0.0 

240.000 

0.0 


2 

0.0 

a 

0 

0 

0 

0 

0 

0 

180.000 

60.CO0 

0.0 


0 

0.0 

10 

0 

0 

0 

0 

0 

0 

180.000 

240. OCC 

0.0 


2 

0.0 

ii 

i 

1 

1 

T 

1 

1 

-100. COO 

240.000 

0.0 


0 

0.0 

12 * 

i 

1 

1 

1 

1 

1 

300. con 

240.00ft 

0.0 


0 

0.0 

GENERATED 

NODAL 

DATA 










NODE 

BOUNDARY 

CONDITION 

CODES 

/-- 

NO CAL POINT 

COORDINATES — 


■/ 



NUMBEP 

X 

Y 

Z 

XX 

YY 

7,7, 

X 

Y 


z 


T 

1 

1 

1 


-1 

-1 

1 

0.0 

0.0 

0.0 



0.0 

2 

1 

1 


-1 

-1 

1 

180.000 

0.0 

0.0 



0.0 

3 

0 

0 


-1 

-1 

0 

0.0 

60.000 

0.0 



0.0 

4 

0 

0 

-1 

-1 

-1 

0 

IRC. 000 

60.000 

0.0 



0.0 

5 

0 

0 


-1 

-1 

0 

0.0 

120.000 

0.0 



0.0 

6 

0 

0 


-1 

-1 

0 

180. 000 

120.000 

0.0 



0.0 

7 

0 

0 


-1 

-1 

0 

0.0 

IflO.OOO 

c.o 



0.0 

0 

0 

0 


-1 

-1 

0 

180.000 

18C.000 

0.0 



0.0 

9 

0 

ft 

- 1 

-1 

-1 

0 

0.0 

240.000 

0.0 



0.0 

10 

0 

6 

-1 

-1 

-1 

0 

180.000 

240.000 

0.0 



0.0 

11 

1 

i 

1 

1 

1 

1 

-10ft. 000 

240.000 

0.0 



0.0 

12 

1 

i 

1 

1 

1 

1 

300.000 

240.000 

0.0 



0.0 

PQUATTON 

NUMBERS 










V 

X 

Y 

Z 

XX 

TY 

zz 







1 

0 

o 

0 

0 

ft 

0 







2 

0 

0 

0 

0 

ft 

0 







.3 

1 

2 

0 

c 

0 

3 







4 

u 

5 

0 

0 

0 

6 







5 

7 

R 

0 

0 

ft 

9 







6 

10 

11 

0 

0 

ft 

12 














Computer Printout 





(Input 

data. 

the 

initial 

design and the 

final design only 

are 

reproduced. ) 


T # W 
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13 

14 

0 

0 

0 

15 

8 

16 

17 

0 

c 

0 

18 

9 

19 

20 

0 

0 

0 

21 

10 

22 

23 

0 

0 

0 

24 

11 

0 

0 

o 

0 

0 

0 

12 

0 

0 

0 

0 

0 

0 


00 
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THREE DIMENSIONAL DEAN ELEMENTS 


NUMBER OF EEAN ELEMENTS = 9 
CONSTRUCTION CODF = 1 
NUMBER OF MATERIALS = 1 
NUNDEP OP GEOMETRIC PROPERTIES- 1 
NUNBEP OF FIXED-END FORCE SETS = 0 


MATERIAL PROPERTY CARDS 


MA T ERT M. 
NUMBER 


SPECIFIC 

WEIGHT 


YOUNGS 

MODULUS 


POISSONS / ALLOWABLE STRESSES — / 

PATIO TENSION COMPRESSION SHEAR 


1 


0. 1C00D 01 0. 3000E OR 0.0 


0.5000E OS 0.S000F OS 0.2885E OS 


GEOMETRIC PROPERTY CARDS 


PROPEFTY 

NUMBER 

X-SECT 

NODE 

X- SECT 

AREA 

X-AXIS Y-AXTS Z-AXIS 

-/ 

1 

1 

0.1000D 01 

0.0 0.0 0.100CE Cl 

0.0 0.0 0.0 

0.0 0.0 0.0 

MOMENTS OF INFPTIA 
SECT MODULI FOP POINT A 
SECT MODULI FOR POINT B 


ELEMENT 

LOAD 

MULTIPLIERS 






A 

B 

c 

D 

X-DTR 

0.0 

0.0 

0.0 

0.0 


T-DIR 

0.0 

0.0 

0.0 

0.0 


Z-DIR 

0.0 

0.0 

0.0 

0.0 



PROCESSED ELEMENT DATA 


ELEMENT 

/—NODE NOS 


/--ELEMENT TD 

NOS-/ 

DESIGN 

VAR 

FIXED 

FND- 

FORCE 

ID 

END RELEASE 

CODES 

BAND. 

NUHBFP 

I 

J 

K 

MATL 

GEONY 

D VAR 

FRACTION 

A 

B 

C 

D 

I 

J 

WIDTH 

1 

1 

3 

12 

1 

1 

1 

0. 1C00E 

01 

C 

0 

0 

0 

cooooo 

oooooo 

3 

2 

2 

U 

12 

1 

1 

1 

0. I^OOE 

01 

0 

0 

0 

0 

000000 

oooooo 

3 

3 

3 

S 

12 

1 

1 

1 

0. 1C00E 

01 

0 

0 

0 

0 

000000 

oooooo 

9 

a 

u 

f> 

12 

1 

1 

1 

O.IPOOE 

01 

0 

0 

0 

0 

000000 

oooooo 

9 

5 

5 

7 

12 

1 

1 

1 

0.1000E 

01 

n 

0 

0 

0 

000000 

oooooo 

9 

fi 

6 

R 

12 

1 

1 

1 

0. 1000E 

01 

0 

0 

c 

0 

oooooo 

oooooo 

9 

7 

7 

9 

12 

1 

1 

1 

0.10C0F 

0 1 

0 

0 

0 

0 

0000^0 

oocooo 

9 

P 

B 

10 

12 

1 

1 

1 

0.100PF 

01 

0 

0 . 

0 

c 

oooooo 

oooooo 

9 

9 

9 

10 

1 

1 

1 

2 

0. 1O00E 

Cl 

0 

0 

0 

c 

oooooo 

0OOOOU 

A 


M.1.9 



BOUNDARY ELEMENTS 


NUMBER OF ELEMENTS = 2 

ELEMENT LOAD MULTIPLIERS 

A B C D 

0.0 0.0 0.0 0.0 

BOUNDARY* ELEMENT DATA 


CONST 

NODE 

/--NODES 

DEFINING 

CONSTRAINT 

DIRECTION--/ 


CODES 

NUMBER 

N 

N I 

NJ 

NK 

NL 

KD 

KR 

1 

q 

11 

0 

0 

0 

1 

0 

2 

10 

12 

0 

0 

0 

1 

0 


STRUCTURE 

STRUCTURE 

LOAD MULTIPLIERS 

LOAD CASE 

A B 

C D 

1 

0.0 0.0 

0.0 0.0 


BUCKLING CONTROL DATA 

COEF»T = 1.00000 

MODEIN = 1 

NMODE = 2 

INDEX = 1 

NVEC = 2 

ALPA = 0.50000 

OMEGA = 0 . 00000 


NODAL POINT LOADS 
NODE LOAD 

NO. CASF RX 

9 1 0.0 

10 1 0.0 


DESIGN VARIABLE INPUT DATA 
DESIGN 

VARIABLE INITIAL MIN ALLOWABLE 

NUMBER VALUE VALUE 

1 0.1500E 01 0.0 

2 0. 1500E 00 0.0 


TOTAL NUMBER OF F0UATTCNS = 24 
BANDWIDTH = 9 
NUMBER OF EQUATIONS TN A BLOCK = 22 
NUMPFR OF BLOCKS = 2 


APPLIED LOADS 

RT PZ MX MT 

0. 100 D 05 0.0 0.0 0.0 

0.1 00D 05 0,0 0.0 0.0 


DISPL ROTATION STTFF 

D R S 

0.0 0.0 U.OOD 01 

0.0 0.0 4.00D 01 




MZ 

0.0 

0.0 


M.'l.lO 



****************************** 
ANALYSIS OF DESIGN NUMBER 0 
****************************** 


NODAL DISPLACEMENTS AND ROTATIONS 


NODE 

NO. 

LOAD 

CASE 

X 

T 


7 

XX 

YT 77 

12 

1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

11 

1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

10 

1 

4.406E-19 

-5. 333E-02 

0.0 

0.0 

0.0 

-5.0772E-21 

o 

1 

4.416E-19 

-5.333E-02 

0.0 

0.0 

0.0 

-1.32SUE-20 

8 

1 

2.004E-19 

-4 • 000E-02 

9.C 

o.o 

0.0 

-3, 017 IE-21 

7 

1 

-7.U67E-20 

-4.9P0E-02 

0.0 

0.0 

0.0 

-4.5553E-21 

6 

1 

6.796E-20 

-2.667E-02 

0.0 

o.n 

o.c 

- 1 . 4P42E-21 

5 

1 

-1.768F-19 

-2 . 667E-02 

0.0 

0.0 

0.0 

5.5310E-22 

4 

1 

1. 172E-20 

- 1 • 333E-02 

0.0 

0.0 

o 

o 

-4.7849E-22 

3 

1 

-8 . 0 10F-20 

-1.333E-02 

0.0 

0.0 

0.0 

2.0715E-21 

2 

1 

0.0 

0.0 

0.0 

o 

o 

0.0 

0.0 

1 

1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

VALUES OP 

DESIGN V APT ABLES 







1 

2 

3 

u 

5 

6 


0 0.1S00F 01 0.1500E 00 


ANALYSIS OF BEAM ELEMENTS , CONSTRN C0DE = 1 


ELEMENT 

X-SECT AREA 

LOAD COND 

AXIAL RX 

SHEAF RY 

SHEA? F7 

TORQUE NX 

MOMENT MT 

MOMENT MZ 

1 

0.1500E 01 

1 

0.1000E 

05 

0.4UP7F-16 

O.C 

o.o 

0.0 

0.2900E-14 




-0. 1000 F 

0 ^ 

-0.44 07F- 1 6 

O.C 

0.0 

0.0 

-0. 2074P-1 5 

2 

0.1S00E 01 

1 

C.1000P 

05 

0.6590E-17 

0.0 

0.0 

0.0 

-0. 1612E-15 




-0. 1000 F 

05 

-0. 559FF- r, 

0.0 

c.o 

0.0 

0.5566F-15 

3 

0.1600F 01 

1 

0. 1000 F 

C5 

0.4487?- 16 

0.0 

o.c 

0.0 

0.207UE-15 




-C.1000E 

05 

-C.4Un7?-if 

0.0 

0.0 

C.O 

0. 2486F-14 

4 

O.I^OOF 01 

1 

0. 1000E 

05 

C. 8590F-17 

0.0 

o.c 

O.C 

-O.ssfif w-15 




-o.nooE 

C5 

-C. 659 OF -17 

0.0 

0.0 

0.0 

0.9520E-15 

5 

0.15C"F 01 

1 

0.1000E 

05 

".4U87E-16 

r >. 0 

0.0 

O.C 

-0. 24R6E-14 




-0. 10C0E 

0 C 

-0.4UR7E-16 

0." 

0.0 

0.0 

0. M77S-14 

6 

0.1500E 01 

1 

C. 1000? 

rs 

C » 6690?- 17 

o.c 

0.0 

O.C 

- 0. 9S20E- 1 5 




-0. 1000c 

o c 

- r .6 CQ C F-17 

0.0 

n. n 

0.0 

0. 7347E-14 

7 

0.1500E 01 

1 

0 . 1 0 00 r 

C 5 

0. 4UR7F- 16 

o.c 

0.0 

0.0 

-0. 517RF-14 
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a 0 . 1500 F 01 


1 



BUCKLING HODF SHAPES 


NOT)F MOPE X T 7 . XX 

NO. SHAPE 


12 

1 

0.0 

0.0 

0.0 

0.0 

0.0 


2 

0.0 

0.0 

0.0 

0.0 

0.0 

11 

1 

0.0 

0.0 

o.n 

0.0 

0.0 


2 

0.0 

0.0 

0.0 

0.0 

0.0 

10 

1 

- 9 . 08 PE -92 

S. 8258-06 

0.0 

0.0 

0.0 


2 

4 . 304 E -05 

- 6 . 652 E -22 

0.0 

0.0 

0.0 

9 

1 

- 9 .CRPE -02 

- 5 . 8258-36 

0.0 

0.0 

0.0 


2 

- 4 . 304 E -05 

1 . 362 E -21 

0.0 

0.0 

0.0 

fl 

1 

- 5 .O 59 F -02 

4 . 36 PE -96 

0.0 

0.0 

0.0 


2 

- 4 , 025 E - 02 

- 4 . 989 E -22 

0.0 

0.0 

o.o 

7 

1 

-S.o 1 GE -02 

- 4 . 36 Q E-D 6 

0.0 

o.c 

0.0 


2 

4 . 025 E -02 

1 . 021 E -2 1 

0.0 

0.0 

0.0 

6 

1 

- 1 . 797 E -02 

2 . 813 E -06 

0.0 

0.0 

0.0 


2 

- 4 .S 31 E -02 

- 3 . 326 E -22 

0.0 

0.0 

0.0 

s 

1 

- 1 . 742 F -02 

- 2 . 913 E -06 

3.0 

0.0 

0.0 


2 

4 . S 3 IF - 02 

6 . 80 PE -22 

0.0 

0.0 

0.0 

4 

1 

- 2 . 3348-03 

1 . 456 E -06 

0.0 

0.0 

0.0 


2 

- 1 . 821 P -02 

- 1 . 663 E -22 

0.0 

0.0 

o.o 

3 

1 

- 2 . 1 14 E -03 

- 1 . 456 E -06 

0.0 

0.0 

0.0 


2 

1 . 82 IE -02 

3 . 404 E -22 

0.0 

0.0 

c.o 

2 

1 

0.0 

0.0 

0.0 

c.o 

0.0 


2 

0.0 

0.0 

0.0 

0.0 

0.0 

1 

1 

0.0 

0.0 

0.0 

0.0 

o.n 


2 

0.0 

0.0 

0.0 

0.0 

0.0 


EVALUATION OF DESIGN NOHSER 


STRESS PATTO LOAD COUP DES VAPIAPJ,F 
MAX 0 . 1333 P GO 1 1 

BIN 9 . 3628 F -20 1 2 


0.0 

0.0 

0 . 7 B 70 E - 14 

0.0 

0.0 

- 0 . 1347 E- 1 U 

0.0 

0.0 

0 . 1743 E -14 

o.o 

0.0 

- 0 . 7 R 62 E- 0 S 

0.0 

0.0 

- 0 . 7862 F -05 


TY ZZ 


o.o 

c.c 

0.0 

0.0 

6 . 5071 E- 04 
-7.8056E-04 
6.6012E-04 
7. 8056E-04 
6.U7R1E-04 
-4. 4184E-04 
6.5293E-04 
4. 4 184E-04 
4 . 1 1 55E-04 
2.6379E-04 
4 .C836E-04 
-2. 637 Q E-04 
1.2060F-04 
5.0896E-04 
1.1446 E-04 
-5.0B96E-04 
0.0 
0.0 
0.0 
0.0 
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PIAX BUCK RATIO? LOAD COND 

0.9S51S 00 1 

0, 6083F 00 1 

DESIGN IS WOT CRITICAL 


DWirOFN SCALING OPEFA^TPN follows 

SCALE FACTOR IS 0.R79ANB OETE n nTNFD BT BHCKLTNG 


CONSTRAINTS 
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****************************** 
ANALYSIS OF DESIGN NTTMBFR U 
****************************** 


NODAL DISPLACEMENTS AND POTATIONS 


NODE 

NO. 

LOAD 
C ASF 

X 

Y 


12 

1 

0.0 

o . n 

3.0 

11 

1 

0.0 

0.0 

0.0 

10 

1 

8.674F- 19 

-9.510E-32 

9.0 

9 

1 

8 . 7 22F- 1 9 

-9.S10E-02 

0.0 

P 

1 

7. 318E-19 

-7. 133P-02 

0.0 

7 

1 

-2.923E-1R 

-7. 133E-92 

3.0 

6 

1 

4. 337F-19 

-4.755E-02 

o 

o 

S 

1 

-2. R17F-1R 

-4.755E-02 

0.0 

4 

1 

1.3S5E-19 

-2.37BE-02 

3.0 

3 

1 

-1 .0R3E-18 

-2. 37BE-02 

3.0 

2 

1 

0." 

0 • 0 

3.0 

1 

1 

0.0 

0.0 

3.0 


7 XX YY 7 Z 


0.0 

0.0 

0.0 

0.0 

0.0 

o 

o 

0.0 

o 

o 

0.0 

0.0 

c 

o 

- 1. 0BU2E- 19 

0.0 

0.0 

-4 . 06 50 E-21 

0.0 

0.0 

-2.4417E-20 

0.0 

o 

o 

-5.U210E-21 

0.0 

o 

o 

2 . 1 6 54 E- 20 

0.0 

0.0 

-4 . 06 50 E- 21 

0.0 

0.0 

2.9793E-20 

0.0 

C.O 

© 

© 

o 

o 

0.0 

C.O 


VALUES OF 

DESIGN VAPI 

ABLES 
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6 

7 

P 

3 0. 

8 4 1 2 E 00 0. 

2R9?p 

no 








ANALYSIS 

OF REAM ELEMENTS , 

CONSTFN C0DE= 

1 






ELEflEN' 1 ’ X 

:-S p CT APEA 

LOAD 

COND 

AXIAL FX 

SHF AR RY 

SMFAB PZ 

TORQUE MX 

MOMENT MY 

MOMENT HZ 

1 

0.R412F 00 

1 


C.1000E 

ns 

0 . 2S59E- 1 S 

0.0 

0.0 

0.0 

0.2051E-13 





-0. 1000F 

os 

-C.2659F-15 

0. 0 

n.n 

0.0 

-0. U554E-14 

2 

0.8412E 00 

1 


0. 1000E 

05 

-C. 1Q00P-16 

0.0 

0.0 

o.c 

-0. 22ROE-14 





-0.1000F 

C5 

0. 1900E-16 

C.O 

o.c 

0.0 

0. 1140E-14 

3 

0.84 1 2E 00 

1 


0.1300E 

05 

0. 2S59E-15 

0.0 

3.0 

0.0 

C. 455UE-14 





-0. 1000 E 

C s 

-C. 265op.-i5 

0.0 

0.0 

0.0 

0. 1140F-13 

4 

0.R412F 00 

1 


n. 1 300P 

05 

-0. 190 Op- IS 

o.c 

0.0 

0.0 

-0. 1 1 40 p - 1 4 





-0.1000F 

fjs 

0. 190CE-16 

O.C 

0.0 

0.0 

-0. 24 4RP-20 

s 

n.R412E 00 

1 


0. 13Q0P 

05 

0. 26SQP- 15 

r. c 

0 . 0 

n .C 

-0. 11U0F-13 





-0.1300E 

05 

-0. 2659^-1^ 

0.0 

O.C 

o.o 

0.2735E-13 

fi 

0.R412E 00 

1 


0.19C0F 

05 

-0. 190CE-16 

o.c 

O.C 

o.c 

-0.743RF‘20 





-0. nocp 

r z 

C. 190OP- is 

o. n 

0 . ° 

0.0 

-0. 114QF-14 

7 

0.8412E 00 

1 


n . 19 n C p 

05 

C. 26*'9F- 1 S 

0. 0 

0.0 

0.0 

-C.2735E-13 
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-0. 1000E OS 

-0. 2659F- 1 5 

0.0 

0.0 

0.0 

0.4331F-13 

0.8412E 00 1 

0.1000E 05 

-0. 1900E-16 

0.0 

0.0 

0.0 

0. 1140E-14 


- C . 1 0 COE 05 

0. 190PE-16 

C.O 

0.0 

0.0 

-0.228CE-14 

0.2R92E 00 1 

0.2310F-15 

-0.2828E-06 

0.0 

0.0 

0.0 

-0.2568F-05 


-0.2310E-15 

0 • 2 R28F- 06 

0.0 

0.0 

0.0 

-0. 2568E-05 


ANALYSTS OF ROU NDAFY FLEMFNTS - CONSTRAINT FOFCFS 
CONST NUMBER LOAD CASE FORCE MOMENT 


1 

1 

-0. 348R6F-16 

0.0 

2 

1 

0. 34694E-16 

0.0 


RUCKLING LOAD PARAMETERS 

0. Q 9R73D 00 0 • 99 95 1 D 00 

RUCKLTNG MODE SHAPES 

NODE MODE X T 

7 

XX 

TY ZZ 

NO. SHARE 
12 1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


2 

0.0 

0.0 

0.0 

0.0 

0,0 

0.0 

11 

1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

10 

1 

-1.563F-C2 

3. 401E-06 

0.0 

0.0 

0.0 

R.223RE-04 


2 

-7.473E-02 

1.827E-05 

0.0 

o.c 

0.0 

3.923PE-04 

9 

1 

-1,56 1 E-02 

-3.401E-06 

0.0 

0.0 

o.c 

-5 . 9975E-04 


2 

-7. 474F-02 

-1 .627E-05 

0.0 

0.0 

0.0 

6.7293E-04 

8 

1 

3. 132E-02 

2. 551E-06 

0.0 

0.0 

0.0 

5.9173E-04 


2 

-4.384P-02 

1. 221E-05 

0.0 

0.0 

0.0 

5.B019E-04 

7 

1 

-4 . fi 44 E-02 

-2.551R-96 

0.0 

0.0 

0.0 

- 3 . 1 1 97E-04 


2 

-2.8S0E-02 

-1.221E-05 

0.0 

o.o 

0.0 

7 , 58 46 E-04 

6 

1 

4.407E-02 

1.700*5-06 

0.0 

0.0 

0.0 

-1.7824E-04 


2 

-1.300E-02 

8. 137E-06 

0.0 

0.0 

0.0 

3.9149E-C4 

5 

1 

-4.5B0E-02 

-1.700 E-06 

0.0 

0.0 

0.0 

3 .21 28E-04 


2 

4.730E-03 

-8. 137E-06 

o.c 

o.c 

0.0 

2.9295R-04 

4 

1 

1 ,«QOr ! *02 

8 . 5Q2E-0 7 

o.c 

0.0 

0.0 

-5.166RE-0U 


2 

1.83 IE-04 

U.C68F-06 

0.0 

0.0 

0.0 

6. 1516E-0S 

3 

1 

-1.7ME-02 

-8.502F-07 

0.0 

0.0 

0.0 

5.C046E-0U 


2 

7.364P-C3 

-4 . 068S-06 

0.0 

o.c 

0.0 

-1. 3913E-04 

2 

1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1 

1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


2 

0.0 

o , 0 

0.0 

O.o 

0.0 

0.0 


******************************* 

EVALUATION OF DESIGN NUMBER U 

******************************* 


STRFSS RATIO 
MAX 0.237PE DO 

MIN D.1B97E-19 


LOAD COND DFS VAPTABLF 
1 1 

1 2 
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NAX BnCK PATIOS 


LOAD COND 


0.1003F 01 1 

0.1000E 01 1 


DESIGN TS CRITICAL 


STRUCTURAL WET GHT = 0.U558E 03 


PEDESTGN OPERATION FOLLOWS 

OPTT PI A LITY INPFX OF DFSTGN VARIABLES FOR BUCKLING CONSTRAINTS 
DV NO ACT/PAS INDEX 

1 ACT 0.100BRE 01 

2 AC T 0.np95RE 00 

NO. OF ACTIVE BUCKLING CONSTRAINTS APE 2 

BUCKLING - CRITICAL DESIGN HAS CONYERGFD 
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N. PLANE STRESS ELEMENTS 


N. 1 Cantilever Beam of Variable Width 



Figure N.1.1 

Cantilever Beam Showing Element and Node Numbers 


The cantilever beam in Fig. N.1.1 is modelled with ten membrane 
elements, the thickness of each element being an independent design 
variable. The results of DESAP 2 for this particular problem can be 
checked, because an exact analytical expression for the width of an 
optimal cantilever beam is available [15]. 

The material properties were taken as 
E = 10 x 10^ psi (Young’s modulus), 
v = 0 (Poisson's ratio), 
a* = a* = 10 6 psi (allowable stress), 
p =1.0 lb./cu. in. (specific weight). 


I 



N.1.2 


The allowable stress was deliberately given a very high value, since 
we did not want the stress constraints to be active. Factor of safety 
against buckling is taken as one, i.e., p* = 1 is used. 

A uniform thickness of 0,1 in, was chosen as the initial design, 
and no lower limits on the thickness were used. The design history is 
summarized in Table N.1.1. The final design, which was reached in 
four iterations, can be seen to compare very well with the exact 
analytical solution, also given in Table N.1.1. 

The solution for an optimal cantilever beam is somewhat difficult 
to extract from Ref. [15]. In the absence of minimum size constraints, 
the optimal thickness distribution can be shown to be 


t(x) 




(N.1.1) 


where P is the desired buckling load, £ is the length of the beam, 
and b represents the fixed depth of the cross section. Inserting the 
properties of our particular problem, the above formula reduces to 


t(x) = 0.24 (1 -x 2 /£ 2 ). 


(N.1.2) 


The thicknesses listed in Table N.1.1 are the values obtained from 
(N.1.2) at the center of each finite element. 

Special notes on input -output : 

1) KSCALE = 1 in Design Control Data informs the program that uniform 
scaling is an exact operation, i.e., the scaled structure does not 
have to be reanalyzed. Note that the stiffness matrix of each 
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element has the form [IL] = [lc^] (t^ = thickness of element) . 

2) The layout of elements lends itself well to automatic generation 

of element and node data compare computer printout of input data 

with the echo of the data cards. 

3) The loading was prescribed on the data card of element No. 10 

see Processed Element Data. Also note the proper use of Element 
Load Fractions (multipliers) and Structural Load Multipliers. 

4) NMODE = 2 in Buckling Control Data asks for analysis of two buckling 
modes , which also means that both of the modes are considered in 
redesign. Since the design is determined by one mode only, the same 
results would be obtained with NMODE =1, but the buckling analysis 
for the second mode would be lost. 

5) INDET = 1 in Buckling Control Data declares the internal forces of 
the prebuckling state to be statically determinate. Consequently, 
the geometric stiffness matrix of the structure does not have to 
be updated after each redesign cycle. 

6) It is important in this problem not to suppress the incompatible 
displacement modes of the elements. The incompatible modes enable 
each element to undergo pure bending deformation, thereby greatly 
increasing the accuracy of buckling analysis, (see Element Control 
Card for plane stress elements in Vol. 1). 

7) The design was terminated when the Optimality Index of each element 
was sufficiently close to one (see Evaluation of Design No. 4) . 



N. 1.4 


Element 

Critical , Scaled Designs 
(thickness in inches) 

Exact 
Solution 
Ref. [15] 

0 

1 

2 

3 

4 

1 

.1939 

.2348 

.2357 

.2381 

.2389 

.2394 

2 

,.1939 

.2283 

.2309 

.2333 

.2341 

.2346 

3 

’.1939 

.2159 

.2213 

.2237 

.2244 

.2250 

4 

.1939 

.1989 

.2067 

.2092 

.2100 

.21Q6 

5 

.1939 

.1789 

.1872 

.1900 

.1908 

.1914 

6 

.1939 

.1578 

.1625 

.1660 

.1668 

.1674 

7 

.1939 

.1378 

.1331 

.1371 

.1380 

.1386 

8 

.1939 

.1207 

.1010 

.1032 

.1043 

.1050 


.1939 

.1084 

.0715 

.0650 

.0659 

.0666 

■S 

.1939 

.1019 

.0531 

.0313 

.0242 

.0234 

Wt (lb) 

96.94 

84.17 

80.14 

79.84 

79.86 

80.00 


Table N.1.1 

Design History of Element Thicknesses 
and Total Structural Weight. 
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CANTTLEVFF B*A M-BI1CKLT W ANALYSTS MB DESIGN USING 11NIPOPM BEMBPANE ELEMENT 


NUMBER ON NODAL POINTS - 22 
NUMBER OF ELEMENT TYPES ■ 1 
NIIHBEP OF LOAD CASES * 1 
NUMBER Of DBS. YAP TABLES ■ 10 


DESIGN CONTROL DA^A 

NCYCl ■ 6 

ESCAIE- 1 

DELTA ■ 0.2500F.-C1 

epstl ■ o.noor oo 

LBDCF « 1 


NODAL POTMT INPUT DATA 


NODE 

BOIINOART 

CONDITION 

CODES 

/--• 

—NODAL POINT 

COORDINATES-- 




NUMBEP 

X 

Y 

7 

XX 

YY 

7.7, 

X 

Y 

z 


T 

1 

1 

1 

-1 

-1 

-1 

-1 

0,0 

0,0 

0.0 

0 

0.0 

2 

1 

1 

0 

0 

0 

0 

0.0 

5 . 00c 

0.0 

0 

0.0 

3 

0 

0 

0 

0 

0 

0 

10.000 

0.0 

0.0 

0 

0.0 

21 

0 

0 

0 

0 

0 

0 

lon.ono 

0 , n 

0.0 

2 

0.0 

U 

n 

0 

0 

0 

0 

0 

10.000 

5.000 

0.0 

0 

0.0 

22 

0 

0 

1 

1 

1 

1 

100, TOO 

5.0PO 

0.0 

2 

0.0 

GENFFATED 

NODAL 

DATA 









NODE 

BOUNDARY 

CONDITION 

CODES 

/— 

•---NODAL POINT 

COORDINATES-- 




number 

X 

Y 

7 

XX 

YY 

7,7 

t 

Y 

z 


T 

1 

1 

1 

-1 

-1 

-1 

-1 

0.0 

0.0 

0,0 


0.0 

2 

1 

1 

-1 

-1 

-1 

-1 

o ,0 

5.000 

0.0 


0.0 

3 

0 

0 

-1 

-1 

-1 

-1 

10.000 

0.0 

0.0 


0.0 

Q 

0 

0 

- 1 

-1 

-1 

-1 

10.000 

5.000 

0.0 


0.0 

5 

0 

0 

-1 

-1 

-1 

-1 

20.000 

0.0 

0.0 


0.0 

ft 

c 

0 

-1 

-1 

-1 

- 1 

20.000 

5 , 0*0 

0.0 


0.0 

7 

0 

0 

-1 

-1 

-1 

-1 

30.000 

O.o 

0.0 


0.0 

P 

0 

0 

- 1 

-1 

-1 

-1 

10.000 

5.000 

0.0 


0.0 

S 

0 

0 

-1 

-1 

-1 

-1 

UO ,000 

0.0 

0,0 


0,0 

10 

0 

0 

-1 

-1 

-1 

•1 

uo.ono 

5.000 

0,0 


0.0 

11 

0 

0 

-1 

-1 

-1 

-1 

50.000 

0.0 

0.0 


0,0 

12 

0 

0 

- 1 

-1 

-1 

-1 

50.000 

5.000 

0.0 


0.0 

n 

0 

0 

-1 

-1 

-1 

-1 

fiO.OOO 

O.o 

0.0 


0.0 

1U 

0 

0 

-1 

-1 

-1 

-1 

fiO.POO 

5, OOP 

0.0 


o.o 

15 

0 

0 

-1 

-1 

-1 

-1 

70. o 00 

0.0 

0.0 


0.0 

16 

0 

0 

- 1 

-1 

-1 

-1 

70,000 

5,ooo 

O.n 


o.c 

17 

n 

0 

-1 

-1 

-1 

-1 

RH.OOQ 

0.0 

o.o 


0,0 

1 R 

0 

0 

-1 

-1 

-1 

-1 

BO. POO 

5,000 

o.o 


o.o 

10 

0 

n 

-1 

-1 

-1 

-1 

qo.mo 

o,o 

0.0 


0.0 

20 

0 

0 

- 1 

- 1 

-1 

-1 

R*.POO 

5.000 

0.0 


0.0 

21 

n 

0 

-1 

-1 

-1 

-1 

100.000 

0.0 

0.0 


0,0 

22 

r 

0 

1 

1 

1 

1 

100.00* 

5.0*0 

0,0 


0.0 


kcwa* ton mumpers Computer Printout 

(Input data, the initial design and the final design only are reproduced) 
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H 

X 

Y 

7 

XX 

YT 

7Z 

1 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

3 

1 

2 

0 

0 

0 

0 

U 

3 

a 

0 

0 

0 

0 

5 

5 

6 

0 

0 

0 

0 

6 

7 

n 

0 

0 

0 

0 

7 

9 

10 

0 

o ■ 

0 

0 

R 

11 

12 

0 

0 

0 

o 

9 

13 

1U 

0 

0 

0 

n 

10 

15 

16 

0 

0 

0 

0 

11 

17 

IP 

0 

0 

0 

0 

12 

19 

20 

0 

n 

n 

0 

13 

21 

72 

0 

0 

0 

0 

1U 

23 

24 

0 

0 

n 

0 

15 

25 

26 

o 

0 

0 

0 

16 

27 

2R 

0 

0 

0 

0 

17 

29 

30 

0 

0 

0 

0 

16 

31 

32 

0 

0 

0 

0 

19 

33 

30 

0 

0 

0 

0 

20 

35 

36 

0 

0 

0 

0 

21 

37 

IP 

0 

0 

0 

0 

22 

39 

on 

0 

0 

0 

0 
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/ ALLOWABLE STRESSES / 

TENSION COPPRESSION SHEAR 

0.10C0E 07 0. 1000 E 07 0.0 


PRN7 BAND 
CODE WDTH 

3 H 

3 8 

3 8 

3 8 

3 8 

3 8 

3 8 

3 8 

3 8 

3 8 



00 


DESIGN VARIABLE INPUT DATA 


PESTGN 

VAFIABLF INITIAL HTN ALLOW ARL F 

NIINBEP VALnR VALUE 


1 0. 1000F OP 0. 0 

2 0.1000F 00 C.O 

3 P.1000P 00 0.0 

U 0.1000E 00 0.0 

5 0 , 10 OOF 00 n.o 

fi O.1000 p 00 0.0 

7 0.1000F 00 0.0 

8 0. 1000E 00 0.0 

9 0.10P0F 00 0.0 

10 O.IOPOE 00 0.0 


TOTAL NnNPER OF EQUATIONS = 90 
BANDWIDTH = R 
NUHRFR OF EQUATIONS TN A BLOCK = 90 
NUHRFR OF BLOCKS = 1 



t-* 

to 



ANALYSIS OP DESIGN N(I«0rp 0 


NODAL DISPLACEMENTS AND ROTATIONS 


NODE 

LOAD 

X 

Y 


7 , 

XX 

YY 

NO • 

CAS* 







27 

1 

-1.000E-01 

-2. 110E-13 

0.0 

o 

o 

0.0 

0.0 

21 

1 

-l.npOE-01 

-2. 1108-13 

0.0 

0.0 

0.0 

0.0 

20 

1 

-9 . 000E-02 

-1.P1UE-13 

0.0 

0.0 

0.0 

0.0 

19 

1 

-9.000P-02 

-1.P14S-13 

0.0 

0.0 

0.0 

0.0 

IP 

1 

-fl.CCOE-02 

-1.520E-13 

0.0 

0.0 

0.0 

0.0 

17 

1 

-8.000E-P2 

-1 . 5205-13 

0.0 

o.p 

0.0 

0.0 

IP 

1 

-7.090P-02 

-1 . 231 E- 13 

0.0 

0.0 

0.0 

0.0 

1 5 

1 

-7 , 000E-Q2 

-1. 23 IE-1 3 

0.0 

O.o 

0.0 

0.0 

14 

1 

-6. 000F-02 

-o. 5*> IE- 1 U 

0.0 

0.0 

0.0 

0.0 

13 

1 

-6.000P-02 

-9.S51F-14 

0.0 

0.0 

0.0 

0.0 

12 

1 

-5 • 0 OOP- 02 

-5.9945-14 

0.0 

0.0 

0.0 

o.O 

11 

1 

-*j« 0 OOP- 02 

-6,9945-14 

0.0 

0.0 

0.0 

0.0 

10 

1 

-4.000P-P2 

-4.7145-14 

0.0 

0.0 

0.0 

0.0 

q 

1 

-4.000E-02 

-4.714F-14 

0.0 

0.0 

0.0 

0.0 

P 

1 

-3.000E-P2 

- 2. 7 86 P - 1 4 . 

0.0 

0.0 

0.0 

0.0 

7 

1 

-3.000E-02 

-2 . 7fl6p- 1 4 

0.0 

0.0 

P .0 

0.0 

6 

1 

-2 ,0005-02 

- 1 . 2965-14 

o.ft 

0.0 

0.0 

0.0 

5 

1 

- 2 . O 00 F - P 2 

- 1 , 2965-14 

0.0 

0.0 

c.o 

0.0 

4 

1 

-1.00PF-02 

-3. 36RP-15 

o 

o 

0.0 

0.0 

0.0 

3 

1 

-1 • 090 E -02 

-3.368E-16 

0.0 

0.0 

0.0 

0.0 

2 

1 

o 

o 

0.0 

0.0 

0.0 

0.0 

0.0 

1 

1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

VALUES OF 

DESIGN V APT AB! ES 







1 

2 

3 

4 

5 

6 


N.1.10 



P 0.10D0F 00 0. 1P03S 00 0. 10P0F 00 0.1Q00F OP O.IOCOE 00 0.1000E 00 O.IOPOF 00 0.1000E 00 0.1000E 00 0.1000E 00 


ANALYSTS OP MFM BP ^ NE FLEETS, C0NST9N CODE= 2 

SIIFET LOAD /---PFAPPAFE FOPCES T N I.OCAI. COORDS //- -H FN BRA N E FORCES IN MATERIAL COORDS-/ 


ELEMENT 

THICKNFSS 

COND 

LOCATION 

N XX 


NY Y 

NX Y 

Nil 


N22 

N12 

1 

0, 1 COPE 

00 

1 

CRN 

-0. 10C0E 

04 

-0.1PP7F-13 

-0.1473E-03 

-0.1C00E 

04 

-0.1007E-13 

-0.1473E-03 

2 

0, 1000F 

00 

1 

CFN 

-O.IOCOE 

nu 

0. 1917F- 14 

-0. 3914F-C3 

-0. 1000F 

04 

0. 1417E-14 

-0.391 4E-03 

3 

0.1P00E 

00 

1 

CEN 

-0. 1PO0E 

04 

0.U968E-13 

-0.4B08F-03 

-0. 1000E 

04 

0.4S6SE-13 

-0.4RORE-03 

u 

o. ipor* 

PC 

1 

CEN 

-0. 1000E 

04 

0.1 156F-12 

-0.3R40E-C3 

-0.10O0F 

04 

0. 1156E-12 

-0.3«40E-C3 

5 

O.IOCOE 

00 

1 

CEN 

-0. 1000F 

04 

0. 1364E-12 

-0.3R91E-02 

-O.IOCOE 

04 

0. 1364E-12 

-0.38B1E-02 

6 

0. 1000F. 

00 

1 

CEN 

-0. 1 P 00 E 

04 

P.2056E-12 

-0.37^4E-02 

-O.IOCOE 

04 

0.2056E-12 

-0.3794E-02 

7 

0.1000E 

00 

1 

CEN 

-0. 10 OPE 

04 

0. 1758F-12 

-0.2419E-02 

-0. 1CC0E 

04 

0. 175RE-12 

-0.2419E-02 

9 

0. 1000* 

00 

1 

CEN 

-0.1000E 

04 

- 0 . Q 5 15F- 1 4 

-O.4O06E-02 

-0. 1000E 

04 

-0.9515E-14 

-0.4006E-02 

Q 

0.1000F 

00 

1 

CEN 

-O.IOCOE 

04 

-0.5750E-13 

-0.2R62E-02 

-O.IOCOE 

04 

-C.5750E-13 

-0.2862E-02 

10 

0 . 10 00 F 

00 

1 

CEN 

0 . 4785 F- 1 3 

-0. 1 00 OE 04 

C.2921E-02 

0.47R5E-13 

-O.IOCOE 04 

0.292 IE-02 


SUCKLING LOAD PAPAHETEB5 


0.51S77D 00 

0.47693D 01 


BUCKLING NODE SHAPFS 







NODF 

rone 

X 

Y 


7. 

XX 

YY 

7Z 

NO. SHAPE 








22 

i 

5.013E-C3 

-1.275E-01 

O.n 

O.p 

0.0 

O.o 



2 

5. 1 32E-03 

4. 278E-02 

7.0 

O.p 

0.0 

0.0 


21 

1 

-5.C13F-03 

-1.275E-01 

0.0 

0.0 

o.o 

0.0 



2 

-5. 132F-03 

4.278E-02 

0.0 

0.0 

o.o 

o.c 


20 

1 

4 , °5 1E-0 3 

-1 • P75E-01 

0.0 

0.0 

0.0 

0.0 



2 

4.570E-03 

6. 22PE-02 

0.0 

0.0 

0.0 

0.0 


19 

1 

-4.9ME-03 

- 1 . 076 F -01 

0.0 

n.o 

n.o 

0.0 



2 

-4.57PE-P3 

6.228E-02 

0.0 

0.0 

0.0 

0.0 


18 

1 

4 .768 F -03 

-9. 807E-02 

0.0 

O.p 

0.0 

0.0 



2 

3.P09E-03 

7.7S1E-02 

0.0 

n.O 

n.o 

0.0 


17 

1 

-4 . 768E- n 3 

-9 . 807E-C2 

0.0 

0.0 

0.0 

0.0 



2 

-3 . 00°E- 03 

7.7S1E-02 

0.0 

o.o 

0.0 

o.o 


16 

1 

4 . 467F-0 3 

-6.9S9F-02 

0.0 

0.0 

0.0 

0.0 



2 

7.R12E-04 

9. 514E-02 

0.0 

o.c 

0.0 

0.0 


15 

1 

-4 . 467E-03 

-6.959E-02 

0.0 

0.0 

0.0 

0.0 



2 

-7.412P-C4 

8. 51UF-02 

0.0 

0.0 

n.o 

0.0 


14 

1 

U,0*6F-03 

-9. 254F-02 

0.0 

0.0 

0.0 

0.0 



2 

-1.69SF-03 

.9,3537-02 

O.n 

0.0 

0.0 

0.0 


13 

1 

-4 , 0 t; 6E-03 

-5.254F-02 

0.0 

o.n 

0.0 

0.0 



2 

1 • 5°5E-03 

R. 3537-02 

o.c 

0.0 

n.o 

n.o 


12 

1 

3.945^-03 

-3 . 73^7-02 

0.0 

o.n 

0.0 

n.o 



2 

-3.629F-03 

7. 303E-02 

o.o 

n.o 

O.n 

0.0 


11 

1 

-3.54RF-03 

-3.713E-02 

0.0 

0.0 

O.n 

0.0 



2 

3. 62 Q F-03 

7 . 3^37-02 

O.p 

0.0 

O.n 

0.0 


10 

1 

2.946F-03 

-2. 4 3 4 E - 0 2 

0.0 

O.p 

0.0 

n.n 



2 

-4 . P6PE-03 

5.595F-02 

O.P 

O.n 

0.0 

0.0 


9 

1 

-2 . 946 E -03 

-2.4347-02 

0.0 

0.0 

0.0 

n.o 



2 

4. P68F-03 

5°5E-92 

0.0 

n , n 

o.o 

n.n 


8 

1 

2. 276F-03 

-1 . 39QF-02 

0.0 

0.0 

n.o 

0.0 



2 

-6 . 046F-03 

3.6037-02 

O.p 

n.n 

n.o 

0.0 
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7 


1 -2 . 2 76 E-03 -1.3R9F-0? 

2 S.Q46E-C3 3.603E-92 

6 1 1.S4QR-03 -S.23RE-93 

2 -4.129E-03 1.760E-92 

5 1 -1.S49E-03 -S.23«E-03 

2 4.12RE-03 1 . 760E-0 2 

4 1 7.R42E-0U -1.56^-01 

2 -2.3 1 UE-03 4 . 6 S 1 E-9 3 

3 1 -7.R42E-04 -I.SfiOE-93 

2 2 . 3 14E-03 4.6S1E-03 

2 1 0.0 0.0 

2 0.0 9.0 

1 1 C.O 0.0 

2 0.^ 7.0 


9.0 

0.0 

0.0 

0.0 

0.0 

o.o 

9.0 

0.0 

o.n 

9.0 

0.0 

o.o 

9.0 

D.O 

0.0 

■9.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

9.0 

rt .o 

0.0 

9.0 

o.o 

0.0 

9.0 

0.0 

0.0 

9.0 

0.0 

O.D 


******************************* 
EVA I.nfcTJ ON OF HESTON NUMBER 0 
***** ************************** 


STRESS PATIO LOAD COND DES VARIABLE 
MAX o . mnOE-M 1 1 

MTN 0.1000P-01 1 10 

MAX BUCK PATIOS LOAD CONE 

0. 1 0 F 01 1 

C.2C07E 00 1 

UNTFOFM SCALING OPEPATION FOLLOWS 

SCAT E FACTOR IS 1.939AND DETERMINED BY BUCKLING CONSTRAINTS 

DESIGN VARIABLES OF SCALED (CRITICAL) DESIGN APF 

VALDFS OF DESIGN VARIABLES 

1 2 3 4 S 

0 0. 1939E CC O.I^qr DO 0. 1939E 00 c* . 19 00 0. 1°39E 00 

STRUCTURAL WFTGHT= 0.969UE 02 

REDESIGN OPERATION FOILOWS 

OPTIMALITY INDEX OF DESIGN VARIABLES FOR BUCKLTNG CONST R A T NTS 
DV NO ACT/PAS INDEX 


1 

ACT 

0. 1 3 7S2E 01 

2 

ACT 

0. 126DRE 01 

3 

ACT 

0. 113R2F 01 


0.0 

0.0 

0.0 

0.0 

0.0 

C.O 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 


6 7 R 9 10 

0.1939E 00 0. 1 9 ?9E 00 0.1939E 09 0.1939E 00 0. 1°39E 00 
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a 

ACT 

0.969S1E On 

5 

ACT 

0.77117? 00 

6 

ACT 

C.*>5 26?E 00 

7 

ACT 

0.364?OE 0 Q 

R 

ACT 

0. 19SS7E CO 

0 

ACT 

0.729<nF-01 

10 

ACT 

0.BS512E-0? 


NO. OF ACTIVE BUCKLING CONSTRAINTS 


AFF 1 


N.1,13 



ANALYSES OF DESIGN NUNpF" 4 

****************************** 


NODAL DISPLACEMENTS AND DOTATIONS 


NODE 

NO. 

Load 

CASE 

X 

Y 


7. 

XX 

22 

1 

-1.0235-01 

-3.9575-14 

0.0 

0 

0 

0.0 

21 

1 

-1.0235-01 

-3.9575-14 

0.0 

".0 

0.0 

20 

1 

-6."93F-02 

-1. 3975-14 

O.o 

0.0 

n.o 

IP 

1 

-6.093F-02 

-3 , 997K- 1 4 

0. 0 

0.0 

O.o 

18 

1 

-4.5755-02 

-2. 8425-14 

0.0 

O.p 

0.0 

17 

1 

-4.575F-02 

-7. 8425-14 

0.0 

0.0 

0.0 

16 

1 

-3.6165-02 

-2. 296F-14 

0.0 

0.0 

0.0 

15 

1 

-3.6 16F-02 

-2. 2965-14 

0.0 

o.o 

0.0 

14 

1 

-2.89 IE-02 

-1 . 7725-14 

0.0 

0.0 

0.0 

13 

1 

-2.8915-02 

-1. 7725-14 

0.0 

O.P 

0.0 

12 

1 

-2.2925-02 

-1. 2995-14 

0.0 

0.0 

0.0 

1 1 

1 

- 2. 2 Q 2E- 0 2 

-1 . 293E-14 

0.0 

0.0 

0.0 

10 

1 

-1 .7685-02 

-8.708F-15 

0.0 

0.0 

O.o 

9 

1 

-1.7685-07 

-9 ,7085-15 

0.0 

9 . 0 

0.0 

P 

1 

-1.291F-02 

-5. 1375-15 

0.0 

0.0 

0.0 

7 

1 

-1 .291E-02 

-5. 13TR-15 

0.0 

0.0 

n.o 

6 

1 

-8. a59E-03 

-2.3P0E-15 

0 .0 

0.0 

0.0 

5 

1 

-B.459F-03 

-2.3905-15 

0.0 

0.0 

0.0 

a 

1 

-4. 187E-03 

-6.7455-16 

o.c 

0.0 

O.o 

3 

1 

-4. 1875-0? 

-6. 246 5-16 

0.0 

0.0 

0.0 

2 

1 

0.0 

0.0 

0.0 

0.0 

P . 0 

1 

1 

0. f) 

n . 0 

0.0 

0.0 

0.0 

VALUES OF 

DESIGN v APT A RL5S 






1 

7 

9 

4 

r. 


YY 


7.Z 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.c 

0.0 

r.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

6 7 a 9 10 
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0 0.23R9E 00 0.2341E OP 0.2244E 00 0.2100F 00 0. 190RE 00 0.1068E 00 0.13ROE 00 0.1043E 00 0.fiS85E-01 0.241RE-01 

ANALYSIS OF MEMBRANE ELEMENTS , CONSTRN COD p = 2 

SHEET T.n\D / MEMBRANE FOPCFS TN LOCAL COORDS //--MEMBRANE FORCES TN MATERIAL COORDS-/ 


ELEMENT 

THICKNESS 

COND 

LOCATION 

NXX 


NYT 

NX Y 

Nil 


N22 

N12 

1 

0 . 2 3 8 Q F 00 

1 

CRN 

-0 . lO^np 

nu 

0 .403SE- 13 

-0. 1602E-03 

-0.1000E 

04 

0.40 35E- 1 3 

-0.1 602E-03 

2 

0.2341E 00 

1 

CEN 

-0. 10COF 

04 

0 . 2325E-13 

-0.2B07E-03 

-0. 1 CODE 

04 

0. 2325E- 1 3 

-O.26O7E-03 

1 

0.2244E 00 

1 

CEN 

-0. 10C0E 

n«i 

-n. ishoe-13 

-0.4477E-03 

-0. 1000E 

04 

-0. 1590E-13 

-0.4477E-03 

4 

0.2100E 00 

1 

CEN 

-C. 1000 E 

C 4 

-0. 1 17 IE- 1 3 

-0.421 IE-03 

-0. 1000E 

04 

-0. 1171E-13 

-0.421 1E-0 3 

5 

0. 1 9 08 E 00 

1 

CEN 

-0. 1009E 

04 

-0.2293F-13 

-0. 3463E-03 

-0. 1 000 E 

04 

-0.2293E-13 

-0.3463E-03 

6 

O.lffRE 00 

1 

CEN 

-0. 1000S 

04 

0.7751E-13 

-0.26U9E-02 

-0. 1000E 

04 

0 • 775 IF- 1 3 

-0 . 2649E-02 

7 

0.138OE 00 

1 

CEN 
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OU 

0 . 1CS Q F- 1 2 
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-0. 1000E 
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8 
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1 
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-0. 1000E 

04 
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0.323BE-14 
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04 

0.3236E-14 

-0.4662E-03 

10 

0.241RE-C1 

1 

CFN 

-0.S914E- 

-14 

-0. 1PO0E 04 

-0.6347E-04 

-0.6914E- 

■14 

-0.1000E 04 

-0.6347E-04 
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22 

1 

-f .097E-03 

1. 227E-01 

o.c 

0.0 

C.O 

0.0 


2 

8. 742E-03 

4 . 7U3F-02 

0.0 

0.0 

0.0 

0.0 

21 

1 

f> .007^-03 

1. 227F-01 

o.c 

o.o 

0.0 

0.0 


2 

-8.74 2E- 03 

4.743E-02 

0.0 

0.0 

0.0 

0.0 

20 

1 

-5.S18E-03 

9.P35E-02 

0.0 

o.o 

0.0 

0.0 


2 

4. S0°?.-03 

7 . 462E- 02 

0.0 

0.0 

C.o 

0.0 

19 

1 

5 . 5 IRE-0 3 

9. 93SE-02 

0.0 

0,0 

0.0 

o.o 


2 

-4. SORE-03 

7 . UB2E- 02 

0.0 

0.0 

0.0 

c.o 

18 

1 

-4.903E-O3 

7. 84RE-02 

0.0 

0.0 

0.0 

0.0 


2 

8. 3S6E-04 

8.S41E-02 

0.0 

0.0 

o.o 

0.0 

17 

1 

4 . 903E-03 

7. R4RE-02 

0.0 

0.0 

0.0 

0.0 


2 

-8.3S6E-C4 

8. 541E-02 

0.0 

0.0 

0.0 

0.0 

If 

1 

-4 . 2R°F-0 3 

6.007F-02 

0.0 

0.0 

0.0 

0.0 


2 

-1.796F-03 

R.34RE-02 

0.0 

0.0 

0.0 

0.0 

IS 

1 

4.289E-03 

6.007E-02 

o.c 

0.0 

0.0 

o.o 


2 

1 .7°fiE-03 

R. 34RE-02 

0.0 

0.0 

0.0 

0.0 

14 

1 

-3 . f 76E-03 

4.413K-02 

o.o 

0.0 

0.0 

0.0 


2 

-3.44QF-03 

7. 294E-02 

0.0 

0.0 

0.0 

o.o 

n 

1 

3.676F-03 

4. 4HE-Q2 

O.P 

0.0 

C.o 

0.0 


2 

3 .449F-03 

7. 294F-Q2 

0.0 

0.0 

0.0 

0.0 

12 

1 

-3. Q63E-03 

3 . Of 4 F-0 ? 

o.o 

0.0 

0.0 

o.o 


2 

-4. 23SE-03 

S.7S1E-02 

0.0 

0.0 

0.0 

0.0 

11 

1 

3. Of 3E-03 

3.064E-02 

o.o 

0.0 

o.o 

o.o 


2 

4.23SF-03 

S . 75 IE- 02 

3.0 

0.0 

0.0 

0.0 

10 

1 

-2.4S0E-03 

1. 961F-D2 

0.0 

c. n 

c.o 

0.0 


2 

-4. 2RSE-03 

4.042F.-02 

0.0 

0.0 

0.0 

0.0 

9 

1 

2. USOd-03 

1 . Q f IE-02 

o.o 

o.o 

0,0 

o.o 


2 

4.2RSF-03 

4. 042F-02 

0.0 

O.o 

0.0 

o.o 

8 

1 

-1 .R3RE-"3 

1. 103B-92 

o.o 

0.0 

o.o 

0.0 


2 

-3. 74CE-03 

2. 432E-92 

0.0 

n .o 

O.o 

0.0 
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1 1.B3RF-0? 1 . 10. IK -02 7.0 0.0 

2 3.740F-03 2 . 4 3 2E -7 2 7.0 O.o 

6 1 - 1 . 225E- 03 4.903E-07 0.0 0.0 

7 -2.7H5F-03 1.171F-02 0.0 0.0 

S 1 1 . 22^E-Q3 4.903E-03 7.0 0.0 

2 2.745E-03 1.1 1 1E-0 2 7.0 0.0 

4 1 -6.126E-P4 1.226E-03 7.0 0.0 

2 -1.440E-03 2.904E-03 7.0 0.0 

3 1 6.126E-04 1.226E-03 0.0 0.0 

2 1.440E-03 2.704E-03 7.0 0.0 

2 1 0.0 0.7 3. n - 0.0 

2 0.0 7.? 7.0 0.0 

1 1 0.0 7 . C 0.0 0.0 

2 0.0 7.0 7.0 0.0 


***** ************************** 
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0 . 100 1 F 01 1 

0. 1ft 1 1E 00 1 
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STRUCTURAL «ET(?HT= 0.79B6E 02 


REDESIGN OPFPATION FOLLOWS 


OPTIMALITY TNOEX OF HESTON VARIABLES FOR BUCK LING CONSTRAINTS 

OV NO ACT/PAS TNHET 

1 ACT o. 10712E 01 

? ACT 0. 177 12E 0 1 

3 ACT 0. 10014F 71 

4 ACT p. 10016F P 1 

5 ACT 0. 1 072 n F 71 

6 ACT o. lOP^fiP PI 

7 ACT 0. 10737P 71 

B ACT 0. 1 70 60 F 71 

9 ACT 0. IP 12 IE 01 

IP ACT P. 0 3 2° 1 F 00 

NO. OF ACTtvf BUCK LI HO CONSTRAINTS AP P 1 

BUCKLT NO - CRITICAL HESTON HAS CONVERGED 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.7 

0.0 

7.0 

c.o 


0.0 

0.0 

0.0 

c.o 

c.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
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0. SHEAR PANEL ELEMENTS 

0 . 1 Cantilever Beam of Pl a te-Girder Construction 

Figure 0.1.1 shows a plate girder modelled with bars and shear 

panels. This problem is quite similar to the cantilever beam in 

Sec. N.l, except that now the shear deformation of the web has a 

significant effect on the buckling load. 

Each shear panel is sized independently, but symmetry conditions 

are imposed on the bar elements that make up the flanges: A^ = A^, 

A« = A., A- n = A„. In addition, all the transverse bars are to 

3 4 19 20 

be of the same size, i.e., A^ = A 22 = ... = A^q. The design variable 
numbers of the elements can be found in Table 0.1.1. 

The design data used in the program is listed below. 


Material properties: 



E(psi) 

V 

a*=a*(psi) 

cr*(psi) 

pClb/in 3 ) 

Shear panels 

200 

0 

- 

3000 

1.0 

Bars 

400 

- 

5000 

- 

1.0 


Element sizes: 



A 

(initial) 

A* 

(minimum) 

Shear panels (in) 

0.1 

0.01 

2 

Flange bars (in ) 

2.0 

0.1 

2 

Transv* bars (in ) 

1.0 

0.1 
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Local buckling of all the elements is to be neglected in the design. 

The design history is summarized in Table 0.1.1. Each design 
was governed by constraints on global buckling and the element sizes, 
i.e., the stress constraints never became active. 

After six redesigns, the convergence criteria are not yet satis- 
fied, although no decrease in structural weight occurred in the last 
three design cycles. The difficulty can be identified by inspection 
of the optimality indices in Evaluation of Design No. 5: the indices 

of all design variables, except that of No. 11 (shear panel No. 1), 
are within the acceptable limits, i.e., only that one element fails 
to satisfy the optimality criterion. 

Additional design cycles would not bring an improvement, because 
we have encountered a fairly common convergence difficulty: oscilla- 

tion of the design about the optimal point due to relaxation factor 
that is too large (a = 0.5 was used). The oscillatory behavior can 
be observed in Table 0.1.1, where design variable No. 11 is seen to 
fluctuate from design No. 3 onwards; an even better indication is 
the optimality index of the element, shown below. 



Design Number 

0 

1 

2 

3 

4 

5 

Opt. Index of 
Des. Var. 11 

0.0301 

0.1517 

0.6261 

1.3781 

0.7076 

1.2672 


The design would converge readily if the procedure recommended in 


Sec. D.l of Vol. 1 is followed. It was suggested that the "normal" 
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value of the relaxation factor a (in our case a = 0.5) is to be used 
for a few design cycles, after which the designs are to be evaluated 
so that appropriate changes to a can be made. In the current problem, 
the oscillatory character of the solution could be diagnosed after 
three redesign cycles, calling for an increase of a (under-relaxation). 
With a = 0.7, the problem would probably converge after a few additional 
design cycles. The optimal design, however would have essentially the 
same weight as the final design in Table 0.1.1. 

Two checks of DESAP 2 can be made using this problem. First, an 
analytical solution for the buckling load of the initial design is 
available on p. 132 of Ref. [13]: 


P 


P 

e 

cr l + nP e / (AG) 


( 0 . 1 . 1 ) 


where P is the critical load without shear deformation of the web 
e 

(Euler buckling load), A and G are the cross-sectional area and the 
shear modulus of the web, respectively, and n represents a shape 
factor (n = 1 for our case). For the problem at hand, we get 
= 1.579 lb. and P ^ = 1.364 lb., which compares favorably with 
P c ^ = 2 x 0.6834 = 1.367 lb. obtained from DESAP 2 (see Buckling 
Load Parameters of Design No. 0). 

The second check is obtained by comparing the final cross- 
sectional areas of the flange elements with the analytical solution 
of Ref. [15], as was done in Sec. N.l. For this specific problem. 
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the optimal flange area can be shown to be 

A(x) = 3.125 (l-x 2 /& 2 ) • (0.1.2) 

Since (0.1.2) does not include the weakening effect of shear deforma- 
tion, the cross-sectional areas obtained from DESAP 2 should be some- 
what larger, which they are indeed, as can be seen in Table 0.1.1. 

Special notes on input -output : 

1) NSCALE = 1 in Design Control Data informs the computer that uniform 
scaling is an exact operation, since the stiffness matrix of each 
element has the form [IC] = [k^]A?, where n = 1. 

2) Extensive use was made of the automatic generation of nodal and ele- 
ment data. 

3) Local buckling of bars was eliminated as a design consideration by 
leaving the moments of inertia blank on the Geometric Property Card 
of the truss elements. Since the computer replaces the blanks by 
10^, the Euler buckling load of each element is too high to govern 
the design. 

4) By setting ISU = 0 (see the Boundary Condition Code in Processed 
Element Data) , local buckling of the shear panels was excluded 
from redesign equations. 

5) NM0DE = 2 in Buckling Control Data requests the use of two buckling 
modes in the redesign process. Because only one mode governs, 

an identical design would be obtained with NM0DE = 1, but all in- 
formation pertaining to the second mode would not be available 
in the computer printout. 
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6) INDET = 1 informs the computer that the prebuckling state is 
statically determinate. Consequently, the geometric stiffness 
matrix of the structure does not have to be updated after each 
redesign. 





tO Element numbers of bar elements 

Figure 0.1.1 

Layout of Plate Girder Showing Node and Element Numbers 
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Des. Var. 
Number 

Element 

Numbers 



Critical 

, Scaled Designs 



Analytical* 

0 

1 

2 

3 

4 

5 

6 


1 

1.2 

2.923 

3.045 

3.273 

3.318 

3.346 

3.345 

3.348 

3.117' 


2 

3,4 

2.923 

2.981 

3.189 

3.235 

3.256 

3.260 

3.259 

3.055 

t N 

3 

5,6 

2.923 

2.858 

3.029 

3.087 

3.105 

3.110 

3.109 

2.914 

£ 

•H 

4 

7,8 

2.923 

2.690 

2.801 

2.877 

2.892 

2.897 

2.896 

2.742 

. 

5 

9,10 

2.923 

2.493 

2.509 

2.604 

2.616 

2.622 

2.621 

2.492 

cr 

t/) 

6 

11,12 

2.923 

2.285 

2.162 

2.267 

2.277 

2.284 

2.284 

2.180 


7 

13,14 

2.923 

2.087 

1.778 

1.864 

1.876 

1.884 

1.884 

1.805 

</) 

Jh 

8 

15,16 

2.923 

1.192 

1.395 

1.396 

1.412 

1.421 

1.422 

1.336 

cd 

PQ 

9 

17,18 

2.923 

1.797 

1.074 

0.888 

0.885 

0.895 

0.897 

0.867 


10 

19,20 

2.923 

1.733 

0.887 

0.493 

0.340 

0.310 

0.309 

0.305 


21 

21-30 

1.461 

0.862 

0.431 

0.216 

0.108 

0.100 

0.100 

- 

/ — s 

11 

1 

1.461 " 

.0888 

.0512 

.0416 

.0495 

.0422 

.0479 

- 


12 

2 

1.461 

.1092 

.0989 

.0935 

.0955 

.0950 

.0950 

- 

•H 

v — / 

13 

3 

1.461 

.1479 

.1545 

.1434 

.1456 

.1451 

.1450 

- 

W) 

14 

4 

1.461 

.2012 

.2065 

.1939 

.1958 

.1951 

.1950 

- 

o 

15 

5 

1.461 

.2639 

.2566 

.2448 

.2459 

.2452 

.2449 

- 

§ 

16 

6 

1.461 

.3298 

.3053 

.2962 

.2960 

.2953 

.2949 

- 

Cu 

17 

7 

1.461 

.3925 

.3507 

.3480 

.3461 

.3454 

.3449 

- 

fn 

cd 

18 

8 

1.461 

.4458 

.3896 

.3991 

.3963 

.3955 

.3949 

- 

O 

19 

9 

1.461 

.4458 

.4185 

.4450 

.4458 

.4457 

.4450 

- 

CO 

20 

1 io 

1.461 

.4845 

.4338 

.4745 

.4854 

.4933 

.4949 

- 

Weight (lb.) 

164.4 

135.5 

121.5 

119.0 

117.8 

117.9 

117.9 

- 


^Neglects shear deformation 


Table 0,1.1 

Design History of Element Sizes and Total Structural Weight 
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BtTCFLING DFSTGN - SHEAR PANEL AND TRUSS ELEMENT 
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NODE 
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NUNBFF 
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XX 

TY 

7.7 
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1 

i 

- 1 

-1 

-1 

-1 

0.0 

-0.500 

0.0 


0 

0.0 

? 

1 


0 

0 

0 

0 

0.0 

0. 5 n 0 

0.0 


n 

0.0 

3 

0 

0 

0 

p 

0 

0 

2 . 500 

-0. 500 

0.0 


0 

0.0 

21 

c 

0 

P 

n 

0 

0 

25. 00° 

-O.50D 

0.0 


2 

0.0 

a 

0 

0 

0 

P 

0 

0 

2.500 

0.500 

0.0 


0 

0.0 

22 

0 

0 

0 

0 

0 

0 

25.000 

0. 5 n 0 

0.0 


2 

0.0 


GENF’i A TED 
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DATA 








NODE 
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C^DES 


NOPAL POINT 
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NUNPSF 

X 

Y 

Z 

XX 

YY 

77 

X 

Y 

7 m 

T 

1 

1 

1 

-1 


-1 

- 1 

0.0 

-0.50C 

0.0 

0.0 

2 

1 

1 

- 1 


-1 

-1 

0.0 

0. 5pn 

0.0 

0.0 

3 

0 

0 

-1 


-1 

-1 

2.50" 

-0.50C 

0.0 

0.0 

a 

0 

0 

- 1 

- 1 

-1 

-1 

2.500 

0. 500 

0.0 

0.0 

5 

0 

o 

-1 

-1 

• 1 

-1 

c .ooo 

-0.5PC 

0.0 

0.0 

6 

0 

n 

-1 


-1 

-1 

5,00^ 

O.5O0 

0.0 

0.0 

7 

0 

0 

-1 

-1 

-1 

-1 

T.^00 

-0.500 

0.0 

0.0 

8 

0 

0 

-1 


-1 

-1 

7.500 

O.sno 

0.0 

9.0 

9 

0 

0 

-1 

-1 

-1 

-1 

10. "0" 

-0.500 

0.0 

0.0 

10 

o 

o. 

-1 


-1 

-1 

10.090 

0.500 

0.0 

0.0 

11 

0 

0 

-1 


-1 

-1 

12.500 

-O.spo 

c.o 

0.0 

12 

0 

0 

- 1 


-1 

- 1 

12.500 

O.soo 

0.0 

0.0 

13 

0 

0 

-1 


-1 

-1 

1 5 .000 

-0.500 

0.0 

0.0 

Itt 

o 

0 

-1 


-1 

-1 

15.000 

0.500 

0.0 

0.0 

15 

0 

0 

-1 

-1 

-1 

-1 

17.500 

-0. 50 0 

0.0 

9.0 

16 

0 

0 

- 1 


- 1 

-1 

17.500 

0.900 

0.0 

0.0 

17 

0 

0 

-1 


-1 

-1 

20.000 

-0.900 

0.0 

0.0 

19 

o 

0 

-1 


-1 

-1 

2C.00O 

0.500 

0.0 

0.0 

19 

0 

0 

-1 


-1 

-1 

22.900 

-o. ^f>0 

0.0 

0.0 

20 

o 

n 

- 1 

- 1 

-1 

- 1 

22 . 5"0 

0. 900 

".0 

0.0 

21 

0 

0 

-1 


-1 

-1 

26.0^0 

-0.5*0 

o. 0 

9 . p 

22 

n 

0 

- 1 

- 1 

-1 

-1 

25.000 

0.9P0 

".0 

0. 0 


EQUATION NUUDFPS 


Computer Printout 


(Input data, the initial design and the last two designs only are reproduced) 
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31 

32 

0 

0 

n 

0 

i q 

3 3 

39 

0 

0 

0 

0 

20 

35 

36 

0 

0 

0 

0 

21 


3R 

0 

0 

c 

0 

22 

39 

90 

0 

0 

0 

p 


0.1.10 



NUMBER OP TRUSS ELEMENTS 


■ 

30 

CONSTRUCTION CODE 


■ 

1 

BOBBER OP WATER TAtS 


* 

1 

BUMPER OP TEMPS POP WHICH PATL 

PROPS 

GIVEN" 

1 

NIJBREP OP DIFFERENT GEOHETRTE3 

PROPS 

GIVEN- 

1 


MATEPIAL 

material 

PROPERTY CARDS 

BHBBPP 5PECTPTC 


YOriNGS 

COEPFT Of /— ALLOWABLE 

STRESSES — / 

NUMBER 

OP TEMPS WEIGHT 

TEMP 

MODULUS 

THFRH PXPAN 

TEMSIOM 

COMPRESSION 

1 

1 0.1*00E 01 

0.0 

O.aonoE 03 

0.0 

0.5000E C4 

0.5000E 04 


GEOMETRIC PROPERTY CAPDS 

GEOMETRY X-5ECT /— MOMENTS OF INERTIA — / 
NUMBER APFA TT 7,7, 

1 0.5000D 00 0.1*0015 07 0, 10001? 07 


ELEMENT 

LOAD 

MULTIPLIERS 

A 

B 

c 

X-DTP 

0.0 

0.0 

0.0 

0.0 

y-dif 

0.0 

0.0 

0.0 

0.0 

Z-DTP 

0.0 

0.0 

0.0 

0.0 

TE**P 

0.0 

0.0 

0.0 

0.0 


PROCESSED BLEB EFT DATA 


ELEMENT 

X 

1 

J€ 

o 

M 

NOS-/ 

/--ELEMENT ID 

NOS-/ 

DESIGN VAF 

REFERENCE 

END FIXITY 

COEFFICIENTS 

RAND 

NUMBER 

I 

J 

NATL 

GEO MY 

n VAR 

FRACTION 

T*NP 

TY 

7.7, 

WID*H 

‘ 1 

1 

3 

1 

1 

1 

0.1000P 01 

0.0 

0.1000D 01 

o.iooon 01 

2 

2 

2 

4 

1 

1 

1 

0.1000E 01 

0.0 

0.10C0D 01 

o.iooon oi 

2 

7 

3 

S 

1 

1 

2 

0.100PP 01 

0,0 

0.1000D 01 

o.iooon oi 

6 

u 

4 

6 

1 

1 

2 

‘ '0. 10 COE 01 

0.0 

0, 10001) 01 

o.iooon 01 

6 

s 

S 

7 

1 

1 

3 

0.1000* 01 

0.0 

o.iooon oi 

0.1000D 01 

6 

6 

6 

A 

1 

1 

3 

0.1C00E 01 

0.0 

0,10000 01 

o.iooon oi 

6 

7 

7 

0 

1 

1 

4 

0. 1000* 01 

0.0 

0.1000P 01 

O.IOOOD 01 

6 

A 

A 

10 

1 

1 

4 

O.IOOOE 01 

0.0 

0.10CCD 01 

O.IOPOD 01 

6 

q 

q 

11 

1 

1 

5 

0.1C0CE 01 

0.0 

0.100PP 01 

o.iooon oi 

6 

10 

10 

12 

1 

1 

S 

o.mro* oi 

0.0 

c.ioocn oi 

O.IOOOD 01 

6 

11 

11 

13 

1 

1 

6 

0.1000* 01 

0.0 

o.iooon oi 

0.10000 oi 

6 

12 

12 

14 

1 

1 

6 

0.1O0PP 01 

0.0 

0.1000D 01 

o.iooon oi 

6 

n 

13 

IS 

1 

1 

7 

0. 10C0E 0 1 

0.0 

0.1000D 01 

o.iooon 01 

6 

iu 

14 

16 

1 

1 

7 

0.10O0P oi 

0.0 

o.iooon oi 

o.iooon oi 

6 

is 

IS 

17 

1 

1 

A 

0.1OCO* 01 

o.o 

o.iooon oi 

o.iooon oi 

6 

16 

16 

Ifl 

1 

1 

8 

O.IOOOE 01 

0.0 

c.ioocn oi 

o.iooon oi 

6 

17 

17 

19 

1 

1 

9 

0.10O0E Cl 

0.0 

o.iooon oi 

o.iooon oi 

6 

1A 

IP 

20 

1 

1 

9 

O.IOOOE 01 

0.0 

C.1000D Cl 

O.IOOOD 01 

f 

iq 

10 

21 

1 

1 

10 

0.10O0F 01 

0.0 

o.iooon ci 

o.iooon 01 

6 

20 

20 

22 

1 

1 

10 

O.IOOOE 01 

0.0 

0.100CD 01 

O.IOOOD 

6 

21 

7 

4 

1 

1 

21 

0.1CP0* 01 

o.o 

o.iooon oi 

o.ircon oi 

4 

27 

S 

6 

1 

1 

21 

O.IOOOE 01 

r .o 

r.ippon oi 

O.IOOOD 01 

4 

23 

7 

a 

1 

1 

21 

/ \inr** ci 

O.o 

o.iooon 01 

o.iooon 71 

a 

74 

0 

10 

1 

1 

21 

n.lPOO* 01 

n .0 

C.ioocn ni 

o.iocon oi 
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0 , 1.11 
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21 
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0. 10COT) 

01 

0.10C0T) 
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u 
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1 

21 

0, 100CE 

01 
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0.1000D 
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0. 1000P 

01 
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1 

1 

21 

0. 10F0F 

01 
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01 
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01 

u 
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21 

0. 100CE 

01 

0.0 

0.10C0D 

01 

0. 1C0QP 

01 

u 
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19 

20 

1 

1 

21 
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0.10C0D 

01 

0. 1000D 

01 

u 

10 

21 

22 

1 

1 

21 

C. 1000* 

01 

0.0 

0.1O00D 

01 

0.10000 

01 

tt 
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HUWPEF OF SHEAR PA NFL ELEMENTS = 10 
CONSTRUCTION COPE * 1 
NUMBER OF MATERIALS = 1 
NUMBER OF TEMPS FOR WHICH NA^L PPOPS GTVFN= 1 


HATFRTAL PPOPERTT CAPPS 




MATERIAL NUMBER SPECIFIC 

NHMPFP OF TEMPS WEIGHT 


TEMP 


YOUNGS 

MODULUS 


n O T SSN 
RATIO 


ALLOWABLE 

SHEAF 


1 


1 0.1000E 01 

o.O 


0.2000E 

03 I 

ELEMENT 

LPAO 

MULTIPLIERS 







A 

B 


c 


X-D T R 

0.0 

O.o 


0.0 


0.0 

T-PTR 

o.O 

0.0 


0.0 


0.0 

7-dir 

0.0 

0.0 


o.O 


0.0 


0.0 


0.3000E 04 


PROCESSED ELEMFNT PA^A 


ELFMFNT 

/ 

— NODE 

NOS 


-EL ID 

NOS-/ 

ROUND 

nES VAR /--EFFECT 

PANEL DIMNS— / 

BAND 

NUMBFR 

T 

J 

K 

L 

(1 ATL 

D VAR 

CODE 

FRACTION 

LONGER 

SHORTEP 

WIDTH 

1 

1 

3 

4 

2 

1 

1 1 

0 

0. 19 OCR 01 

0.0 

0.0 

4 

2 

3 

5 

ft 

4 

1 

12 

0 

0.1000F 01 

0.0 

0.0 

3 

3 

5 

7 

8 

6 

1 

13 

0 

0.1000E 01 

o.c 

0.0 

8 

a 

7 

9 

n 

8 

1 

14 

0 

O.IPPOE 01 

0.0 

0.0 

8 

ft 

Q 

11 

12 

10 

1 

15 

r 

0. 1000* 01 

0.0 

0.0 

3 

ft 

11 

13 

14 

12 

1 

1ft 

n 

0, 1000E 01 

0.0 

0.0 

3 

7 

13 

15 

1ft 

14 

1 

17 

0 

0.10P0E 01 

0.0 

0.0 

8 

8 

15 

17 

13 

1ft 

1 

13 

0 

0.1000E 01 

0.0 

0.0 

8 

Q 

17 

19 

20 

18 

1 

19 

n 

0.1000E 01 

0.0 

0.0 

8 

10 

19 

21 

72 

20 

1 

20 

0 

0. 10C0F 01 

0.0 

0.0 

8 


STRUCTUPE 
LOAD CASE 


STRUCTURE LOAD MULTIPLIERS 


0.0 


0.0 


0.0 


0.0 


BUCKLING CONTROL DATA 


coefft 

NODEIN 

NMOPE 

TNDET 

NVEC 

ALPA 

OMEGA 


1.00000 

0 

2 

1 

? 

0.50000 

0.80000 


NODAL POINT LOADS 


NODE IPAD 
NO. CASE 


”1 


ALLIED LPADS 
RY *Z 
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21 1 - 0 . 1 00 n 'M 0.0 0.0 0.0 0.0 

22 1 - 0 . 100 D M 0,0 0.0 0.0 0.0 


DESIGN V API ABLE TNPTTT DA r A 


DESIGN 




VAPI ABLF 

INITIAL 

HI N ALLOWABLE 

NUtIDFF 

VALUE 

VAL"F 

1 

0.2000E 

01 

0. 1000E 00 

2 

0. 2000 F 

01 

0. 1 OO^F 00 

3 

O.2000F 

01 

0. 1000E 00 

4 

0 . 20 OOF 

01 

0. 1 0 00 E 00 

5 

0.2000F 

01 

C.1000F 00 

5 

0.2000E 

01 

0. 1009F 00 

7 

0 , 20 OOF 

01 

0. 1009E OD 

R 

0 * 20 OOF 

01 

0. 1000E 00 

q 

0.2000F 

01 

0.1000F 00 

10 

0.2QC0F 

01 

0. 1000F 00 

ii 

0.1000E 

00 

0.1000E-01 

12 

0 * 1 0 COE 

00 

0. 1000F-0 1 

13 

0.10C0E 

00 

0. 1°005-Q1 

14 

0. 1QD0F 

00 

0. 10005-0 1 

15 

0.1000F 

CO 

0. 10005-01 

15 

0. 1 0 COE 

00 

0. 1000F-01 

17 

0. 1000F 

00 

C, 1000E-01 

1 R 

0.1000F 

00 

0. 10005-01 

19 

0 . 1 00 CE 

CO 

0. 10005-01 

20 

0.10CCF 

00 

0, 100OR-01 

21 

0. 1 0 OOF 

01 

0.1000E 00 


TOTAL NMHRFP OF EQTIATTONS = 40 
BANDWIDTH = R 
NUHFFP OF EQUATIONS IN A BLOCK * 40 
WONPFP of FLOCKS = 1 


0.0 

0.0 
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ANALYSTS OF DESIGN WINDED 0 


NODAL DT 

SPLACFSFNTS AND NOTATIONS 





NODE LOAD X 

Y 


7 

XX 

TT 

NO. CASE 






22 

1 -3,1258-02 

-1.693E-15 

0,0 

0.0 

0.0 

0,0 

21 

1 -3.1258-02 

-1,6938-15 

0,0 

0.0 

0.0 

0.0 

20 

1 -2.B12F-02 

-1.919E-15 

0.0 

0.0 

0.0 

0.0 

19 

1 -2.B12F-02 

-1,9108-19 

0.0 

0.0 

0.0 

0.0 

IB 

1 -2. 500E- 02 

-1.196F-15 

0,0 

0,0 

0.0 

0.0 

17 

1 -2.600E-02 

-1,1968-15 

0,0 

0.0 

0.0 

0.0 

16 

1 -2,1ft7*-02 

-B.B09K-16 

0,0 

ft, 9 

0.0 

0.0 

IS 

1 -2.1B7R-02 

-ft , R90R-16 

0.0 

o.ft 

0.0 

0,0 

19 

1 -1.R75F-02 

-9,5998-16 

0.0 

O.ft 

0.0 

0.0 

13 

1 -1.875F-02 

-6,5958-16 

0.0 

0.0 

0.0 

0.0 

12 

1 -1.562F-02 

-9,6958-16 

0.0 

0,0 

0.0 

0.0 

11 

1 -1.563R-02 

-9 , 6B5B-16 

0.0 

0.0 

0.0 

0.0 

10 

1 -1.2908-02 

-3. 0908-16 

0.0 

0.0 

0,0 

0.0 

9 

1 -1.250R-02 

-3.090F-16 

0.0 

ft,n 

0.0 

0.0 

ft 

1 -9.3758-03 

-1 , 73**8-16 

0,0 

0.0 

0.0 

0,0 

7 

1 -9.37SE-03 

-1.717E-16 

0,0 

0.0 

0.0 

0.0 

6 

1 -6.250F-03 

-7,7678-17 

0.0 

0.0 

0.0 

0.0 

5 

1 -6.250E-03 

-7,7671-17 

0.0 

0.0 

0.0 

0.0 

9 

1 -3. 125F-03 

-1,9928-17 

0.0 

0.0 

c.o 

0.0 

3 

1 -3.1258-03 

-1,9928-17 

0,0 

0.0 

0.0 

0.0 

2 

1 0,0 

0,0 

0.0 

0.0 

0.0 

c.o 

1 

1 o,o 

0,0 

0.0 

o.ft 

0.0 

o.c 

tALHFS OF DESIGN YADT ADMS 






1 

2 

.1 

9 


6 



tt 


7 A 9 10 
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0.2P00E 
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o. 20 op n 

01 

0. 20 0">F 

01 

0, 2000E 

01 

0.2000F 

01 

0. 20 COE 

01 

0.2000E 

01 

10 

20 

0. 1000F 
0. 1000E 

or 

01 

0. 1009? 

00 

0.1000F 

00 

0 . 1 00 OF 

no 

0.1O00E 

00 

0. 1000 F 

00 

0.1000E 

on 


ANALYSTS 

OF TFIJSS 

ELEMENTS , 

, CONST FN COOE= 1 

ELEMENT 

X-SFCT AREA LOAD 

CON 0 

AT TAT. F^PCF 

1 

0,?OO0E 

01 

1 

-o. moop 01 

2 

0.2C90E 

01 

1 

-C.1000E 01 

3 

0.2000F 

01 

1 

- n . 1090? 01 

M 

0.2000E 

01 

1 

-0.10P0F 01 

i; 

0.2000? 

01 

1 

-0. 1000 E ni 

6 

0. 2000P 

01 

1 

-0. 1 n 00F 01 

7 

0.20COE 

01 

1 

-0, 100C? 01 

R 

0 . 20 OOF 

01 

1 

-0. iron? n 1 

q 

C.2000E 

01 

1 

-0. I^OO* r '^ 

10 

0. 200PF 

01 

1 

-0.10 n 0* ni 

11 

0.2000F 

0 1 

1 

- 0 . inno? n i 

12 

0*2000f 

01 

1 

-0.1000? 01 

1 3 

0 . 20 OOF 

C 1 

1 

-0. 10 00? 01 

in 

0.200"E 

01 

1 

-0.10C9E 01 

15 

0.2000F 

0 1 

1 

-0. Ipnof 01 

16 

9.2P00E 

01 

1 

-0 . 1000 F 1 

17 

0.2000E 

01 

1 

-o. icon? oi 

IP 

0.2000E 

01 

1 

-o. mooE oi 

19 

9.2C90F 

01 

1 

-0. 1000? 01 

20 

0.2090* 

Cl 

1 

-0.1000E 01 

21 

0. imoF 

01 

1 

0. 065OE-1B 

22 

0. 1000E 

01 

1 

-0 . 25 80 ?- 1 7 

23 

0. 100 Op 

01 

1 

-0. 64 1 ? R- 1 7 

24 

0. 100CE 

Cl 

1 

-0. 1BP0E-17 

25 

0. lOnOF 

01 

1 

-n. 9R7f ?- 1 7 

26 

0. 100 OF 

01 

1 

-0. 393CH-17 

27 

0. 1000E 

01 

1 

-0. 1115E-16 

2B 

0. 1 0 OOF 

01 

1 

-n # 19 Q 6 F- 17 

29 

0. 10F0F 

0 1 

1 

-0. 544RS-17 

30 

0. 1000E 

01 

1 

-0. 1077F- 16 


ANALYSTS 

OP SHFAP PANEL? 

CCNSTRN CnOE= 1 










ELEMENT 

THICKNESS 

COND 

t a 

K L 

SHFAF FLOW 

1 

0. 1000? 

00 

1 

-0. 4657F-«.JP -0. 4657r.QR 

-0, 4657F-0P -C.4657E-0R 

-0.4657E-0R 

2 

0. 1000E 

00 

1 

-n, 12HE-07 -0. 1211F-07 

-0. 121 lr-07 -0. 1211F-07 

-0.1211^-07 

3 

0. 10POF 

00 

1 

-0.782 3 F-07 -0.7R23F-0 7 

-F.7R23F-07 -C.7R23F-07 

-0.7R23E-07 

u 

0. 10P0E 

on 

1 

-0.7637F-07 -0.76 3 7 F-07 

-0.7637F-07 -0.7637F-07 

- 0 . 7f 37E -O 7 

5 

0. lOnOF 

on 

1 

-o. *215F-07 -0. e 2l5F-O7 

-0.5215F-07 -C.5215E-C7 

-0.5215E-07 

6 

0. 10CO? 

00 

1 

-0.1910F-06 -0. 1° 1 9F-°6 

-0 . 19 1 Q F-0 6 -0. iqi4F-06 

-0.1919E-06 

7 

0. 1000? 

on 

1 

-0.14030-06 -0. 1900F-06 

- C. 14C0P-06 - 0 . 1 Q OPF - 06 

-0. 1 9 00F-C6 

R 

0. 1F00E 

00 

1 

-0. 1RR1E-06 -0. 1 P R 1 F-06 

-6. 1RP1F-06 -0. 1PP1E-06 

-0.1RB1F-H6 

q 

0. 10 OOF 

00 

1 

-0.1PF3P-06 -0. 106 3 F-06 

-0. 1P6 3F-06 -C.1863F-06 

-0. 1R63E-06 

10 

0. 10 O n P 

on 

1 

-0.1R4UF-06 -0.1R44F-06 

-n.inunr-06 -O.iruuf-CB 

-0.1P44E-06 


0.2909E 

0.1000E 


01 

no 


0. 2000 F 01 O.2C0^F. 01 

0.1000E 00 P.1000? QO 


BUCKT.T NG LOAF papAMFTFRr 


0 . 1.16 



0.68432P 00 


9.29B94n 01 


MJCFT.ING EODF SHAPES 


HOPE HOPE X T 7. v T TT 

RO.SHRPF 


?2 

1 

-ft . 65BE-02 

3. IflfiF 00 

0.0 

O.n 

n.O 

0.0 


2 

4.25^-02 

1.042E 00 

0.0 

0.0 

0.0 

0.0 

21 

1 

ft . 65RE- 02 

3.186E 00 

0.0 

0.0 

o.O 

0.0 


2 

-4.255F-0? 

1.0U2E 00 

O.o 

o.o 

0.0 

0.0 

20 

1 

-ft • 55 IE- 02 

2.6ft«F 00 

0.0 

0.0 

0.0 

0.0 


2 

3.772E-07 

1.552E 00 

o.n 

0.0 

0.0 

0.0 

19 

1 

8.551F-C2 

2.5RPE 00 

0.0 

o.n 

0.0 

0.0 


2 

-3 . 772E-02 

1 . *'52E 0" 

0.0 

n.n 

0.0 

0.0 

1 9 

1 

-ft , 234E- 02 

2.202E on 

0.0 

O.n 

0.0 

0.0 


2 

2. 437E-02 

1.946E 00 

0.0 

O.o 

0.0 

0.0 

17 

1 

R. 2345-02 

2.202E 00 

0.0 

o.O 

0.0 

0.0 


2 

-2.432F-02 

1.9U6P 00 

0.0 

0.0 

0.0 

0.0 


1 

-7.7 1 4E- 02 

1.740E 00 

0.0 

0.0 

o.o 

0.0 


2 

5.714P-03 

2 . 1 1*15 00 

0.0 

O.o 

o.n 

0.0 

15 

1 

7.714F-C2 

1.7U0E 00 

0.0 

0.0 

0.0 

0.0 


2 

-5.9 14 E-0 3 

2. 116E On 

0.0 

0.0 

0.0 

0.0 

14 

1 

-7 . ft 05E- 02 

1.313E 00 

0.0 

0.0 

0.0 

0.0 


2 

-1 . 359E-02 

2.04fiE 00 

0.0 

o.o 

0.0 

0.0 

13 

1 

7 . 005E-02 

1.313E 00 

0.0 

0.0 

0.0 

n.O 


2 

1. 359F-02 

2.04PE no 

0.0 

0.0 

o.n 

0.0 

12 

1 

-5.122F-02 

9. 333E-01 

0.0 

0.0 

0.0 

O.C 


2 

-2.9R5E-02 

1.792E 00 

0.0 

0.0 

0.0 

0.0 

11 

1 

fi. 122E-02 

0 . 333E-01 

0.0 

0.0 

0.0 

0.0 


7 

2.9R5E-02 

1.392E 00 

0.0 

O.n 

0.0 

0.0 

10 

1 

-5.OR9F-02 

5.0 BSE- 01 

0.0 

O.O 

0.0 

0.0 


2 

-3.^547-02 

1.3UUF 00 

0.0 

n.O 

0.0 

0.0 

q 

1 

5.0R9F-02 

5 . 085E- 01 

0.0 

0.0 

0.0 

0.0 


2 

3 • 954F-02 

1. 344E 00 

0.0 

O.n 

0.0 

0.0 

ft 

1 

-3.931E-02 

3.473E-01 

0.0 

0.0 

0.0 

0.0 


2 

-4 . 052E-02 

q. ^ qop-0 1 

0.0 

0.0 

0.0 

0.0 

7 

1 

3.931E-02 

3.473E-01 

0.0 

0.0 

0.0 

o.o 


2 

4,0^2F-02 

ft . 599E-01 

0.0 

0.0 

n.O 

0.0 

6 

1 

-2.fi75E-02 

1 . 5605-01 

o.o 

0.0 

o.n 

0.0 


2 

-3.296E-C2 

4.ng^E-01 

0.0 

O.n 

0.0 

0.0 

5 

1 

2. 635f- 02 

1.5R0E-01 

0.0 

0.0 

O.n 

o.o 


2 

3.2R6F-02 

4. 0R5F-01 

o-.n 

0.0 

O.n 

0.0 

4 

1 

- 1 i 354F- 02 

1.923F-02 

0.0 

0.0 

0.0 

c.o 


2 

-1 .ft27« , -02 

1 . 044 E-0 1 

0.0 

0.0 

O.o 

O.o 

3 

1 

1 • 354**- 02 

1. 923E-02 

0.0 

0.0 

c.o 

c.o 


2 

1.R27F-02 

1.044E-01 

0.0 

0.0 

O.n 

0.0 

2 

1 

0.0 

o.n 

0.0 

n.O 

0.0 

0.0 


2 

0.0 

O.n 

0.0 

0.0 

O.n 

O.C 

1 

1 

0.0 

0.0 

0.0 

O.C 

0.0 

c . i ) 


2 

0.0 

0 . 0 

0.0 

n.O 

0.0 

n.O 


******************************* 

.EVMIIATTOM OF OFSTON NUMBER 0 

******************************* 

STPESS »ATTO LOAD conn PES vaptmilf 


I 



0 . 1.17 


SAX 

NIK 


o.iooor oo 

0.6000F-01 


11 

1 


0 
0 

WAX BUCK PATIOS LOAD CO*D 

0.1461F 01 1 

0.3346* 00 1 


HKTFOnn SCALING OPERATT OK FOLLOWS 


SCALP FACTOR TS 1.461AND nETERSINFP FT BUCKLING CONSTRAINTS 


PESIGN VAPTABLES OF SCALED (CRITICAL) DESIGN APE 


VALUES OF DESIGN VARIABLES 


1 2 


3 4 


5 6 


S 9 10 


0 

10 

20 


0.2923F 01 
0 . 1461F 00 
0 • 146 IE 01 


0.2923S 01 
0.1461* 00 


0.2923F 01 
0 • 146 IE 00 


0.2923E 01 
0.1461E 00 


0.2923E 01 0.2Q23E 01 0.2923E 01 
0, 146115 00 0. 1 461 E 00 0.1461F 00 


0.2923E 01 
0* 1461 F 00 


0.2923E 01 0.2923E 01 

0* 1461 F 00 0.1461E 00 


STRUCTURAL WEIGHT- 0. 1644E 03 


REDESIGN OPERATION FOLLOWS 


OPTTN ALTTT INDEX OF DESIGN VARIABLES FOR BUCKLING CONSTRAINTS 
DV NO ACT/PAS INDEX 


1 

ACT 

0.76543E 00 

2 

ACT 

0.72B20F 00 

3 

ACT 

0.6S71BE 00 

4 

ACT 

C.SS991E 00 

6 

ACT 

0.44533E 00 

6 

ACT 

0.324R4F 00 

7 

ACT 

0* 2 1026 F 00 

9 

ACT 

C.11279E 00 

Q 

ACT 

0.41971F-01 

10 

ACT 

0.47413^-02 

11 

ACT 

0. 10O71E-01 

12 

AC* 

0.26621F 00 

13 

ACT 

0. 01639* 00 

14 

ACT 

• 0.13336F 01 

IS 

AC* 

0.20S04E 01 

16 

ACT 

0.2924SF 01 

17 

ACT 

0. 3*613* 01 

IB 

ACT 

C.4169SE 01 

19 

ACT 

0.461S7f 01 

20 

AC* 

0.49549* 01 

21 

ACT 

■o.o 


NO. OF ACTIVE BUCKLING C7N STR A T NTS APE 1 


0 . 1.18 



****************************+* 

ANALYSIS OF DESIGN NUMBER 5 

** ********* * + + * * * * *********** * 


NODAL DISPLACEMENT S AND ROTATIONS 


NODE 

NO. 

IOAD 

CASE 

X 

Y 


z 

XX 

YT 

22 

1 

-4 .7955*02 

-2. 2435-15 

0.0 

0.0 

0.0 

0.0 

21 

1 

-4.795E-02 

-2.2445-15 

0.0 

o.c 

c.o 

0.0 

20 

1 

-2.7775-02 

-1 .949E-16 

0.0 

0.0 

0.0 

0.0 

19 

1 

-2.7775-02 

-1.9505-15 

0.0 

0.0 

o.o 

o 

o 

18 

1 

-2 . 079E-02 

-1.653E-15 

0.0 

o.c 

o 

o 

0.0 

17 

1 

-2.C7QE-02 

-1.6535-15 

0.0 

o.c 

o.o 

o.o 

16 

1 

-1.639F-02 

-1.366E-15 

0.0 

0.0 

0.0 

o 

o 

IS 

1 

-1. 6392-02 

-1.3675-15 

0,0 

0.0 

0.0 

0.0 

14 

1 

-1.9075-02 

-1. 1025-15 

0.0 

0.0 

0.0 

0.0 

13 

1 

-1.3075-02 

-1.102E-15 

0.0 

0.0 

0.0 

o 

o 

12 

1 

-1 . 03417-02 

-8. 6245-16 

0.0 

0.0 

0.0 

0.0 

11 

1 

-1.0345-02 

-8. 621E-16 

0.0 

0.0 

0.0 

0.0 

10 

1 

-7.9525-03 

-6. 5025-16 

0.0 

0.0 

0.0 

0.0 

9 

1 

-7,9525-03 

-6.4955-16 

0.0 

''.C 

0 . 0 

c.o 

8 

1 

-5.79SE-P3 

-4. 6295-16 

0.0 

0.0 

o.o 

0.0 

7 

1 

-6. 79^5-03 

-4. 629E-16 

0.0 

O.o 

0.0 

0.0 

6 

1 

-3.98SE-D3 

-3 . 07 IE-1 6 

0.0 

0.0 

0,0 

0.0 

5 

1 

-3 . 7855-0 3 

-3.0325-16 

0.0 

0.0 

0.0 

O.o 

4 

1 

-1 .869F-03 

-1 .7185- 16 

0.0 

0.0 

0,0 

0.0 

3 

1 

- 1 . 88 85-03 

-1.71^5-16 

0 .0 

O.C 

0,0 

0.0 

2 

1 

0.0 

0.? 

0.0 

0.0 

c.o 

o.c 

1 

1 

0 ." 

o.o 

0.0 

0.0 

O.o 

0.0 

VALnns of 

DESIGN VARIABLES 







1 

2 

3 

4 

c 

6 


0.1.19 



0 

0.334SE 01 

0. 3260 E 01 

0.31 10E 

Cl 

0. 20*17? 0 1 

0. 2622F 

01 

0.2284E 

01 

0. 1B84E 

01 

9.1421F 

10 

20 

0.4222P-01 
0.100CF 00 

0.9^01 *-01 

0. 14? 1 F 

00 

0. 19 ? IK 00 

0.2452K 

00 

0.2953F 

00 

0.3454? 

00 

P,395?F 


ANALYSTS 

OF Tpnss FLFMFNTS , CONSTFN CODF = 1 


ELEMENT 

X-SECT APFA 

LCAT) COND AXIAL FORCE 


1 

0.3345F 01 

1 

-C.100CF Pi 


2 

0.3345E 01 

1 

-O.IOOOE 01 


3 

0.3260P 01 

1 

-0.1 POOF 01 


4 

0 . 3260 E 01 

1 

-O.IOOOE 01 


5 

0.3110F 01 

1 

-C.1090F 01 


6 

0.3110*=* 01 

1 

-O.IOOOE 01 


7 

0.28*17? 01 

1 

-O.IOCOR 01 


a 

0.2897* 01 

1 

-O.IOOOE 01 


9 

0.2622F 01 

1 

-0.1 POOF 01 


10 

9.2622F 01 

1 

-O.IOOOE 01 


1 1 

0.2284? 01 

1 

-O.IOOOE 01 


12 

0.22R4F 01 

1 

-O.IOOOE 01 


1 3 

0. IPBUp 01 

1 

-O.IOOOE 01 


14 

0. 1994K 01 

1 

-O.IOOOE 01 


1? 

0.1421K 01 

1 

-O.IOOOE 01 


ir> 

0.1421K 01 

1 

-O.IOOOE 01 


17 

0.89R2F 00 

1 

-O.IOOOE 01 


IB 

0.8952F 00 

1 

-O.IOOOE 01 


1° 

0 . 30 97*=- 00 

1 

-O.IOCOF Cl 


20 

0.3097? 00 

1 

-O.IOOOE 01 


21 

o. icon* oo 

1 

0. 134RE-17 


22 

O.IOOOE 00 

1 

0.4 055E- 1 7 


23 

0 . 1 000 F oo 

1 

-0. 1950E- 17 


24 

0 . 1 o 00 K 00 

1 

-0.292 IE-16 


?* 

0.1000F 00 

1 

-0.1 33UE-16 


26 

0.1000E 00 

1 

0.6324F-18 


27 

0 . moor 00 

1 

O.fiP^qF-TT 


2P 

0.1000K oo 

1 

0 • 2 524F- 1 8 


2* 

0 . 1 POPE 00 

1 

0.3862E-16 


30 

O.10O0K 0Q 

1 

0.3166E-16 


an at. rain 

0* SHEAP PANFLS, CONS TR N CODE^ 1 






AVFPAdE 



L J ft 11 


ELEWFNT 

THICKNESS 

COND 

I .1 K L 

SHEAP FLOW 

1 

9.4222E-C1 

1 

-9.531RE-CR - P . 5 3 1 PF-0 R -9.5318F-0R -0.531RE-C8 

-0.531 RE- 08 

2 

0. °50 1*-0 1 

1 

-0.1133E-07 -0. 1133F-07 -0.1133E-07 -0. 1 1 33E-07 

-0. 11 33E-07 

3 

0. 14 51* CO 

1 

-0.7261E-07 -0.7261F-07 -n.726lF-07 -0.7261E-C7 

-0.7261F-07 

4 

0. 1 9 5 1 E 00 

1 

-0.9990E-07 -0 . 9° 90 F- n 7 -O.Q990F-O7 -0.9990F-07 

- 0 . Q 9 90F- 07 

s 

0.2U52E 00 

1 

- 0 . 223RE- 96 -0 . 223P^-06 -0. 223RE-06 -0.223RF-06 

-0. 223RE-06 

6 

0. 2953K 00 

1 

-0.1P31F-06 - n . 1P31F-06 -0.1831^-06 -0. 1B31F-06 

-r. 1 R 3 IF- 06 

7 

0.T4S4E CC 

1 

-0 . 204 1 F- 06 -0 . 2 n 4 1F-06 -0.2941E-06 -0.2041P-06 

-0. 20 41 H- 06 

a 

0.3QS5F CO 

1 

-0.2332F-06 -0.2332F-06 -9.2332F-06 -P.2332E-06 

-0.2332F-06 

Q 

0.4457? 00 

1 

-0. 1366F-06 -O. 1366F-0^ -0. 1366E-05 -0. 1366F-05 

-0. 1366E-05 

1 r 

0.4933F 00 

1 

-O. 19^7F-C? 19?7f-0*> -0.19S7E-05 -C. 1957E-05 

-0. 1957F-05 


0.8952? OP 
0.U457F 00 


0.3097E 00 
0.4933F 00 


niTCKLTNn TOAD PAPAHETFPS 


0 . 1.20 



0.99984n 00 


0. 20937f) 01 


BTTCKr TNG MODE SHADES 


NODE 

no.: 

"ODE 

SHAPE 

X 

Y 


7, 

XT 

22 

1 

1 • 1S6E-C1 

-3. 124E 00 

9.0 

9.0 

C.O 


2 

1.01SE-02 

1.037R 00 

9.0 

o.n 

0.0 

21 

1 

-1 . 1 56E-01 

-3. 1 24 E 00 

9.0 

0.0 

0.0 


2 

-1.01SE-02 

1.037E 00 

o.n 

n.o 

0.0 

2C 

1 

1.0U1F-01 

-2. 5S2F 00 

9.0 

0,0 

0.0 


2 

7.9P1E-03 

1.0B6E 00 

9.0 

n.O 

0.0 

1 9 

1 

-1.0U IE-01 

-2.SS2E 00 

9.0 

n.O 

0.0 


2 

-7.981E-03 

1.CB6F 00 

0.0 

o.n 

0.0 

IB 

1 

9. 248F-02 

-2.037E 00 

0.0 

o.n 

0.0 


2 

5 • Q 0 1E-0 3 

1. 12SE 00 

9.0 

0.0 

0.0 

17 

1 

-9.2U8E-02 

-2.037F 00 

9.0 

0.0 

0.0 


2 

-5. 90 IE-03 

1. 125R 00 

0.0 

O.n 

0.0 

16 

1 

R.991P-02 

-1. PB1E CO 

9.0 

O.C 

0.0 


2 

3.969E-03 

1.1S2K 00 

0.0 

0.0 

0.0 

IS 

1 

-B . 09 1 p-02 

-1.5B1E 00 

9.0 

0.0 

0.0 


2 

-3.969E-03 

1.1S2E 00 

O.C 

0.0 

0.0 

14 

1 

6 • 93SE-02 

-1. 1B2E 00 

9.0 

n.O 

0.0 


2 

2. 197E-03 

1. 1 70E 00 

9.0 

O.n 

0.0 

13 

1 

-6.93SE-02 

-1 . 1B2E 00 

9.0 

0.0 

o.n 


2 

-2. 197E-03 

1 . 170F 00 

0.0 

0.0 

0.0 

12 

1 

5.779E-02 

-R. 4 195-01 

0.0 

0.0 

0.0 


2 

5. 947F-C4 

1. 17BF 00 

9.0 

0.0 

0.0 

11 

1 

-5.779E-02 

-8. U10E-01 

9.0 

0.0 

n .0 


2 

-5. 9U7F-0U 

1 . 17RE 00 

9.0 

0.0 

0.0 

10 

1 

U.623E-02 

-S. S78E-01 

9.0 

n.o 

0.0 


2 

-B.270E-OU 

1. 177E 00 

3.0 

0.0 

0.0 

9 

1 

-4.623E-02 

-S.57BE-01 

9.0 

0.0 

0.0 


2 

B.270F-0U 

1. 177E no 

0.0 

0.0 

0.0 

B 

1 

3.46 Bf-02 

-1. 324E-01 

n.O 

0.0 

0.0 


2 

-2.054E-03 

1. 167P! 00 

9.0 

0.0 

0.0 

7 

1 

-3.46BE-0? 

-3.124P-01 

9.0 

0.0 

0.0 


2 

2.054E-C3 

1 . 16715 nr) 

0.0 

0.0 

o.n 

6 

1 

2. 312F-02 

-1.64RE-01 

O.C 

0.0 

0.0 


2 

- 3. 061 E- 03 

1. 1 40E 09 

9.0 

O.C 

0.0 

S 

1 

-2. 312F-02 

-1. 64BE-01 

9.0 

0.0 

n.o 


2 

3. 06 IE-03 

1 . 14RE 00 

0.0 

0.0 

n.o 

4 

1 

1. 1S7F-02 

-S. 49SE-02 

9.0 

0.0 

0.0 


2 

- 3 . 792E- 03 

1. 115E 90 

9.0 

0.0 

0.0 

3 

1 

-1.1S7F-02 

-S. 49SE-02 

9.0 

0.0 

0.0 


2 

3.792F-C3 

1.115E 00 

9.0 

n.o 

0.0 

2 

1 

0.0 

0.9 

9.n 

n.O 

0.0 


2 

0.0 

0. n 

9.0 

O.n 

O.n 

1 

1 

0.0 

9.0 

9.0 

O.C 

P.0 


2 

0.0 

n.O 

9.0 

O.n 

0.0 


EVALUATION OF DESIGN NH.MRSP 5 

*4 ************** ******* ** ****** 


STRESS P A T T 0 


LOAD OOND DFS V A P I A P T 



0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

n.o 

0.0 

0.0 

0.0 

0.0 

0.0 

n.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0 . 1.21 


MAX 

0.100P p 01 


0 

HTN 

0. 2P27F.-01 


0 


MAX RUCK PATIOS 

LOAD 

COND 


0. 1000F 01 

1 



0 • 41 76 E 00 

1 


DESIGN 

IS CRITICAL 




STFUCTI1PA1 WETGHT= 0.1179* 93 
REDESIGN OPERATION FOLLOWS 

''PTT M A LIT Y T NP EX’ OF DESIGN VARIABLES FOR BUCKLING CONSTRAINTS 
PV NO ACT/PAS TNPEX 


1 

ACT 

0. 10016E 

01 

? 

ACT 

0. 90931 p 

00 

3 

ACT 

0. 90931E 

00 

U 

ACT 

C. 

00 

5 

ACT 

0. 9 09 s OF 

00 

f 

ACT 

0. n qq7q 7. 

00 

7 

ACT 

C. 1D002E 

01 

n 

ACT 

0, 1001 IE 

0 1 

0 

ACT 

0. 10035F 

01 

10 

ACT 

0. nq732E 

00 

ii 

ACT 

0. 12B72P 

01 

12 

ACT 

0. 1000 U E 

01 

13 

ACT 

0. 09P67F 

no 

1U 

ACT 

0. 9DR72F 

00 

15 

ACT 

0.^9796* 

no 

16 

ACT 

0 . n 9 7 fi q 

00 

17 

ACT 

0 . 09753F 

00 

IP 

ACT 

0. 9 9 70 Q P 

on 

in 

ACT 

0.99675E 

00 

20 

AC T 

0 . innfifE 

0i 

21 

PASS 

0.0 



NO. OF ACTTVF BUCKLING CONSTRAINTS ARE 1 


0 . 1.22 



***+***********+************** 

ANALYSTS OF DESIGN HOPPES fi 

* ft**************************** 


NODAL DT5PL ACENFNTS AND POTATIONS 


NODF 

NO. 

LOAD 

CASE 

X 

T 


7 

XX 

TY 

22 

1 

-4. 796E-C2 

-1.777E-1S 

0.0 

0.0 

0.0 

0.0 

21 

1 

-4.796E-0? 

-1.777E-15 

0.0 

0.0 

0.0 

0.0 

20 

1 

-2.776E-02 

-1 . 506E-15 

0.0 

0.0 

0.0 

0.0 

ID 

1 

-2 . 776 E- 02 

-1. 507E-1S 

0 . c 

0.0 

0.0 

0.0 

1 fl 

1 

-2.07QF-O2 

-1.237E-15 

0.0 

0.0 

0.0 

0.0 

17 

1 

-2.079F-02 

- 1 . 2 3RE- 1 5 

0.0 

0.0 

0.0 

0.0 

16 

1 

-1.630^-02 

-9 . RB1E-16 

3.C 

n.n 

0.0 

0.0 

15 

1 

-1.639F-02 

-9.R77F-16 

0.0 

0.0 

0.0 

0.0 

1U 

1 

-1.307E-02 

-7. 594 E- 1 6 

o 

o 

0.0 

o.o 

0.0 

13 

1 

-1 .307E-02 

-7 . 590S- 1 6 

0.0 

0.0 

0.0 

0.0 

12 

1 

-1.034E-0? 

-D.620F-16 

0.0 

0.0 

0.0 

0.0 

11 

1 

- 1 . 034E- 02 

-5.61 RE- 1 6 

0.0 

O.o 

0.0 

0.0 

10 

1 

-7.H53E-03 

-1.9R9E-16 

0." 

0.0 

0.0 

o.c 

9 

1 

-7. 9S.3E-03 

-3. 9R5E-16 

0.0 

0.0 

o.o 

0.0 

P 

1 

- c .7R5E-03 

-2. 671E-16 

0.0 

0.0 

0.0 

0.0 

7 

1 

-5.795E-03 

-2. 669P-16 

o.n 

o.c 

0.0 

0.0 

6 

1 

-3.705F-O3 

-1.619P-16 

0.0 

o.o 

0.0 

0.0 

5 

1 

-3.785P-03 

-1.619E-16 

0.0 

0.0 

C.o 

0.0 

u 

1 

- 1 . B67E-0 3 

-7.R70E-17 

0.0 

0.0 

0.0 

0.0 

3 

1 

-1 . 867 E-03 

-7.R74R-17 

0.0 

o.n 

n.o 

o.c 

2 

1 

0.0 

D.O 

o.n 

n .P 

o.o 

c.o 

1 

1 

0.0 

0,0 

0.0 

o.o 

c.o 

0.0 

VALUES OF 

DESIGN VARIABLES 







1 

2 

3 

4 

Cl 

f. 



.23 


0 0.33Q8E 01 0. 325°7 01 0.31O9F 01 C.2R°6F 01 0.2621E 01 0.2284F 01 0.1884” 01 0.1422E 01 0.8967F 00 0.3093F 00 

10 0.47867-01 P.°503E-01 0.1480E CO O.I^OF HO 0.2449E 00 0.2949F 00 0.3449E 00 0.3949P 00 P.4450F 00 0.494«F, DO 

20 0.1000E 00 


analysts of tpijss ftehfnts, constrn code= 1 

P1.FHFMT T-SECT AB FA LOAD r^ND AXTAL P9PCF 

1 0.3148F 01 

2 0.3348F 01 

3 o.32SOE 01 

4 0.325OE 01 

5 0.31C9F 01 

6 0.3 10PE ni 

7 0 . 2 fi 96F 01 

8 0.2996” 01 

8 0. 26 2 IF 01 

10 n .2621” 01 

11 0.22AUF 01 

12 0.2284F 01 

13 0.1884F. 01 

14 0. 1984F C 1 

16 0. 1422E 01 

16 0. 1U22F 0 1 

17 0.R967E 90 

IP 0.B967E 00 

19 0.30Q3E 00 

20 0.3093F O0 

21 0. 1000F 00 

22 O.ICOOF 00 

23 O.IOOOF 00 

24 O.IOOOF 00 

2* O.ICOOF 00 

26 0.1090E 00 

27 0. 1 C OOF 00 

2P 0.1000E 00 

29 0.1^00? CO 

3 C O.IOOOF 00 


ANALYSTS OF SHF AR PANFLS, CONST BN CODF= 1 

LOAD / SH FA” FLOW AT MODES / AVERAGE 

ELEUPBT THICKNESS COHD I 0 K L SHEAR FLOW 

1 D.U7R6F-01 1 -9 . c 38 9” -OR -0.53PPF-08 -0.538BE-08 -0.S3BRE-0R -C.S38RF-0R 

2 0. 9503F-01 1 - 9 ■ 1304F-07 -0.13O4E-07 -0.130U”-C7 -0.1304F-07 -0.13047-07 

3 0.1460F CO 1 -0.68P3F-C7 -9. 6PR0”-n7 -O.fiPPnE-97 -C.68PCF-07 -0.6RR0F-07 

4 0. 1 q 50F 00 1 -n.RB74F-07 -Q.PR74F-^7 -0.8874F-07 -C.8R74F-07 -0.RP747-07 

6 0.24U9E DO 1 -0 . 224 2F- n 6 -0.2242F-06 -0.224 ?p-06 -0.2242E-06 -0.22427-06 

6 0. 2949E 00 1 -0.2181 E-06 -D.2151F-96 -O.21MF-06 -0.2151F-06 -0.21517-06 

7 0.34UQF no 1 -0. 22?8”-06 -0.772SF-06 -0^22287-06 -0.2228F.-Ofi -0.222SE-06 

R 0.394PE CO 1 -0. 20677-06 -0 . 2067F-06 -0.2967^-06 - C . '>*677-06 -0.2067E-06 

9 p.UUSOF 00 1 -o . IP 1 1 p-05 -0.1911P-08 -0.191 IE-05 -C.IPIIF-O 6 -0119117-0^ 

10 0. 4 94 Pp on 1 -0, 13 Q lF-08 -C.1991F-08 -0.1301”-05 -C.1791F-05 -0.1391E-08 


nnCKLTNG LOAD PARAMETERS 


I 


1 - 0 . 1 onQE 0 1 
1 -0.1O00E 01 
1 -0. 1 000 E 01 
1 -0. 1000F 01 
1 -0.1000E 01 

1 -0. 1000E 01 

1 -0. 1CC0E 01 
1 -0. 1000E 01 
1 -O.innoF ”1 
1 -O.inOOF 01 
1 -0*1000” 01 
1 -0. in 00 E 01 
1 -O.IOOOF. 01 
1 -0. 10COE 01 
1 -O.IOOOF 91 
1 -n. lone” 01 

1 -O.IOOOE 01 
1 -0. 1C00E 01 
1 -O.IOOOF 01 
1 -C. innc-r n 1 
1 0. 1 Q C7 17 
1 -0.9 3 PRE-IP 
1 -0 . fl 1 P 4 F- 1 7 
1 -0. 1692F-16 
1 -0.8269E-17 
1 -0.16997-16 
1 -0. 14057-16 
1 C.87S3E-16 
1 0.2R04F-16 
1 -0.10SOE-1R 


0.1.24 



0.99997p CO 


0 . 236°0D 01 


BflCFLTNfi MODE SHAPES 


NODE BODE 
NO. SHAPE 

X 

Y 


7. 

XT 

22 

1 

1.1 56E- 01 

-3.11 «E 00 

0.0 

0.0 

0.0 


2 

-i.om -02 

-1.061P 00 

o.n 

0.0 

0.0 

21 

1 

- 1 . 1 5f>E- 01 

-3.119F 00 

O.o 

O.n 

O.n 


2 

1.niQE- n 2 

-1.061F 00 

0.0 

0.0 

0.0 

20 

1 

1 . Oft IE-01 

-2.547K 00 

0.0 

0.0 

0.0 


2 

-7.R57F-03 

-1. 1 10H 00 

0.0 

0.0 

0.0 

19 

1 

- 1 . 04 1E-0 1 

-2.547* 00 

0.0 

0.0 

0.0 


2 

7, R57P-03 

-1. 1 10 F OO 

0.0 

n.o 

0.0 

ie 

1 

9.250E-02 

-2.032E 00 

0.0 

n.p 

n.o 


2 

-5.637F-03 

-1. 1URF Of) 

0.0 

0.0 

0.0 

17 

1 

-° . 250F- 02 

-2.032B 00 

0.0 

0.0 

0.0 


2 

5.R77E-03 

-1. 1U8F 00 

0.0 

O.n 

o.n 

Ifi 

1 

R.0O3E-02 

-1.57fiF 00 

0.0 

o.n 

n.o 


2 

-3. SB IE-03 

-1. 1 74 F 00 

0.0 

O.n 

0.0 

15 

1 

-B.093F-02 

-1.57RG 00 

0.0 

0.0 

n.o 


2 

3 . *'R1E-03 

-1 . 174 E 00 

0.0 

0.0 

o.o 

1ft 

1 

fi.937E-02 

-1.177F 00 

0.0 

0.0 

0.0 


2 

-1.705E-03 

-1. 189P 00 

0.0 

0.0 

0.0 

13 

1 

-6.937E-02 

-1. 177E no 

0.0 

0.0 

0.0 


2 

1 . 705E-03 

-1. 1R9E OO 

o.o 

0.0 

o.o 

12 

1 

5.7R1E-02 

-R.357E-01 

0.0 

o.n 

0.0 


2 

-1.931F-OS 

-1. inSE 00 

0.0 

n.o 

0.0 

11 

1 

-5.7P1F-02 

-R. 357E-01 

0.0 

o.n 

o.o 


2 

1.931E-05 

-1.195F 00 

o.n 

0.0 

0.0 

10 

1 

ft, fi2ftF- 02 

-5.524K-01 

o.n 

0.0 

o.o 


2 

1.UR1F-03 

-1. 190K 00 

0.0 

0.0 

0.0 

9 

1 

-4 .62ft p - 02 

-5, 524E-01 

0.0 

o.o 

0.0 


2 

-1 , 4f* 1 E -0 3 

-1 . 1 90 K 00 

o.n 

0.0 

0.0 

R 

1 

3.46PF-02 

-3.270E-01 

0.0 

0.0 

0.0 


2 

2.719F-03 

-1, 177F 00 

0.0 

0.0 

0.0 

7 

1 

-3.46PE-02 

-3.270E-01 

0.0 

0.0 

0.0 


2 

-2.71QF-03 

-1. 177E 00 

0.0 

O.n 

0.0 

ft 

1 

2.312E-02 

-1. 594E-01 

0.0 

0.0 

o.o 


2 

3. 726F-03 

-1. ISftE 00 

0.0 

O.n 

0.0 

5 

1 

-2. 312E-02 

-1.594E-01 

0.0 

0.0 

0.0 


2 

-3.726E-03 

-1. 154F 00. 

0.0 

0.0 

0.0 

ft 

V 

1, 1SSF-02 

-4.OR1E-02 

0.0 

O.n 

0.0 


2 

U.419F-03 

-1. 1 15F 00 

0,0 

0.0 

0.0 

3 

1 

-1.1SSF-02 

- 4 . q fi IE- 02 

0.0 

n.o 

0.0 


2 

-4.419E-03 

-1 . 1 1 s F 00 

0.0 

n.o 

o.n 

2 

1 

0.0 

0.0 

0.0 

n.n 

0.0 


? 

0.0 

n.O 

o.c 

0.0 

0.0 

1 

1 

0.0 

0.0 

0.0 

O.n 

o.n 


2 

0.0 

0.0 

0.0 

0.0 

n.o 
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P.1.1 


P. PLATE ELEMENTS 


P. 1 Circular Plate with Simple Support 



Figure P . 1 . 1 

Simply Supported Plate 
under Uniform Compression 


The problem under consider- 
ation is to find the optimal 
distribution of thickness for 
the plate shown in Fig. P.1.1. 

As the design will be axi- 
symmetric, it is sufficient to 
model only a sector of the plate 
by finite elements, as indi- 
cated by the shaded region in 
the figure. 


The finite element model of the sector is shown in Fig. P.1.3. 

The conditions of symmetry on the x-axis are imposed by the boundary 
condition codes in the nodal point input data, but two sets of boundary 
elements are required on the "skewed" boundary of the segment. 

Boundary element Nos. 1-10 are rotational springs that suppress 
rotations about the boundary line, whereas element Nos. 11-20 re- 
present extensional springs that prevent the normal, in-plane displace- 
ments. The spring constants of the boundary elements must be much 
larger than the corresponding stiffnesses of the elements, illustrated 
in Fig. P.1.2, in order to be effective in suppressing the boundary 


displacements and rotations. We used the ratios ^ound^^plate = 5000 


P.1.2 



lP(ib.) 

U 


,_ fs 



Figure P.1.2 

Calculation of (a) Rotational, and (b) Extensional Stiffness of 
a Plate Element 


for the rotational springs and 500 for the extensional springs, 

approximately, where kp late is computed for the initial design. 

The material properties used in the design are: 

E = 10^ psi (Young's modulus), 

V = 0.3 (Poisson's ratio), 

a* = a* = 50,000 psi (allowable stress), 

3 

p = 1.0 lb. /in (specific weight), 
p* = 1.0 (lower bound on the critical load parameter). 

The design was started with a uniform thickness of 0.5554 in., 
and all of the ten plate elements were sized independently. No 
constraint was placed on the minimum thickness. The design history 
of the plate is given in Table P.1.1. The optimal design, which is 
governed by the buckling constraint alone 9 was reached after five 
redesign cycles. 




P.1.3 


An analytical solution for the optimal plate with a contin- 
uously varying thickness has been derived in Ref. [16]. The results 
are plotted in Fig. P.1.4 together with the discrete element thick- 
nesses obtained from DESAP 2. The thickness distributions correlate 
quite well; the small discrepancy in the two weights can be attributed 
to the continuous vs. discrete variation of thickness. 

Special notes on input -output : 

1) Buckling deformation of a flat plate does not involve stretching 

of the middle surface, i.e., only the bending stiffness matrix of 
each element: [IC] = [k^]A n , n = 3, offers resistance to buckling. 

Consequently, uniform scaling is an exact operation with KSCALE = 3 
(see Design Control Data). 

2) The variation in the spring constants of boundary elements shown 
in Boundary Element Data (higher near the center than at the edge) 
reflects similar changes in the stiffnesses of the plate elements. 
The idea here is to maintain an approximately constant ratio of 

^bound^plate t ^ rou S^ ou ' t the finite element model. 

3) NMODE = 2 in Buckling Control Data requests that two buckling 
modes be considered in the design. Since only one mode governs 
the final design, the results would be identical if NMODE = 1 
were used, but the buckling load and the mode shape of the second 
mode would not appear in the printout. 

4) INDET = 1 in Buckling Control Data informs the computer that the 
internal forces of the prebuckling state are statically determin- 
ate, i.e., the geometric stiffness matrix of the structure does 
not have to be recomputed after each design. 


P.1.4 


5) The slight variation of the membrane forces from uniform, biaxial 
compression (see Analysis of Plate/Shell Elements of Design No. 0) 
is caused by the presence of the boundary elements. 



Plate element 
D Boundary element 



Figure P.1.3 

Finite Element Model for a Two Degree Sector of Plate Showing Node and Element Numbers 
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Element 

Critical (scaled) Designs (thickness in inches) 

0 

1 

2 

3 

4 

5 

i 


.6641 

.7110 

.7218 

.7226 

.7222 

2 

1 

.6540 

.6997 

.7112 

.7121 

.7117 

3 

.5554 

.6345 

.6799 


.6966 

.6969 

4 

.5554 

.6078 

.6513 

.6677 

.6719 

.6731 

5 

.5554 

.5765 

.6137 

.6322 

.6384 

.6406 

6 

.5554 

.5436 

.5671 

.5865 

.5948 

.5982 

7 

.5554 

.5121 

.5123 

.5289 

.5389 

.5436 

8 

.5554 

.4846 

.4535 

.4566 

.4668 

.4729 

9 

.5554 

.4630 

.3995 


.3706 

.3772 

10 

.5554 

.4486 

.3626 


.2530 

.2371 

Wt (lb)* 

96.92 

89.94 

86.68 

85.14 

84.54 

84.59 


*The weight of a 2° sector 


Table P.1.1 


Design History of Element Thicknesses and Structural Weight 





Plate thickness (in. 



Figure P.1.4 

Comparison of the Results of DESAP 2 with the Analytical Solution 
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NUHRFP OF NODAL POINTS = 31 
NUMPFP OF FLFMFNT TYPFS = 2 
NUMRFR O p LOAD C A S c 5 = l 
MIJMPFP. OF OF S . VARtAPL r S = 10 
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DESIGN CONTROL DATA 

NCYCL = 10 

KSC AL F- 3 
DFLTA = 0.2000F-01 

FPSIL = 0.1000F 00 

L PUCK * 1 


NOPAL POINT INPUT DATA 


nodf poumdapy condition rnoES / urcM point rnr°niNATj-<; / 


NUMRFP 

X 

Y 

z 

XX 

YY 

11 

X 

Y 

Z 


T 

1 

1 

-1 

0 

-1 

1 

-1 

n.O 

n.n 

0.0 

0 

0.0 

2 

0 

0 

0 

0 

0 

0 

10.000 

o.o 

0.0 

0 

n.n 

20 

0 

1 

1 

1 

0 

0 

lon.onn 

o.o 

0.0 

2 

o.o 

3 

0 

0 

0 

0 

0 

0 

9.994 

0.349 

0.0 

0 

0.0 

21 

0 

0 

1 

0 

0 

1 

99.940 

3.490 

O.o 

2 

0.0 

22 

1 

1 

1 

1 

1 

1 

A. *04 

inn.?R9 

0.0 

0 

0.0 

31 

1 

l 

1 

1 

1 

1 

*6 .<*50 

103.4*0 

0.0 

1 

n.o 


gfnfratfd nopal data 


NHDF BOUNDAPY CONDITION CODFS / NOTAl POINT r. nrc nf A AT e 5 / 


NIJMBFP 

X 

Y 

7 

XX 

YY 

71 

X 

Y 

l 

T 

1 

1 

-l 

0 

-1 

1 

-1 

n.O 

O.o 

o.O 

o.o 

2 

0 

-1 

0 

-I 

0 

-1 

lo.o no 

0.0 

0.0 

o.n 

3 

0 

0 

0 

0 

0 


9.994 

0. *40 

0.0 

n.o 

4 

0 

-1 

0 

-1 

0 

-1 

20.000 

0.0 

0.0 

0.0 

6 

0 

0 

0 

0 

0 

"l 

19.9R3 

0.498 

0.0 

n.n 

6 

0 

-l 

n 

-\ 

0 

-l 

TC.onn 

n.n 

0.0 

o.n 

7 

o 

0 

0 

n 

0 

-l 

79 .SR2 

1.047 

O.n 

0.0 

fl 

0 

-l 

0 

-1 

0 

-1 

40.000 

n.n 

n.n 

0.0 

9 

0 

0 

0 

0 

0 

“1 

39.976 

l • *9 6 

0.0 

n.o 

10 

0 

-1 

0 

-1 

0 

-1 

so. non 

o.n 

o.o 

n.o 

11 

0 

0 

0 

0 

0 

-1 

49.970 

1.745 

o.n 

0.0 

12 

0 

-1 

n 

-1 

0 

-1 

60.0 00 

0.0 

O.o 

0.0 

13 

0 

0 

0 

0 

0 

-1 

59.9(4 

7.094 

o.o 

o.n 

14 

0 

“1 

0 

-1 

0 

-1 

7C .pnn 

o.o 

o.o 

o.n 

15 

0 

0 

0 

0 

0 

-1 

69.950 

2.44 3 

0.0 

O.n 

I 6 

0 

-1 

0 

-1 

0 

-1 

on. coo 

n.n 

n.n 

n.n 

17 

0 

0 

n 

c 

0 

-1 

79 .p*P 

7. 79 7 

0.0 

O.n 

IB 

0 

-1 

0 

-1 

0 

-1 

90.000 

n.o 

n.o 

n.n 

19 

0 

0 

0 

0 

0 

-1 

PS .94* 

*.14 1 

n.o 

n.C 

20 

0 

1 

1 

1 

0 

-1 

inn. pm 

0.0 

n.n 
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n.n 

n.o 
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0.0 

0.0 

0.0 

0.0 

0.0 

o.n 

n.O 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.n 

0.0 

0.0 

c.o 

0.0 

O.n 

n.n 

0.0 

0.0 

c.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

n.o 

0.0 

0.0 

O.o 

0.0 

0.0 

n.o 

0.0 

1.1616F-03 

-3.3272F- 

0.0 

-1 .0485F-03 

?.00*4F- 

0.0 

0.0 

-3.? ' a 07F- 

0.0 

0.0 

■3 .O066 F - 


P.l.16 




2 

0.0 

0.0 

-2.B68E- 

-01 

l .7003r-03 

-4.R*Q6F-02 

0.0 

e 

1 

0.0 

0.0 

-1.824= 

00 

0.0 

-2.18ROF-P2 

0.0 


? 

0.0 

0.0 

-2.R6RR- 

-01 

0.0 

-4.0 7 26F-02 

0.0 

7 

1 

0.0 

0.0 

-2.019= 

00 

5. 946 2F-04 

- 1 . 703 2F-02 

0.0 


? 

0.0 

0.0 

-7 • 8ftB= - 

-01 

1.731 ?F -03 

-4.0585=-02 

n.O 

6 

1 

0.0 

0.0 

-2.019F 

00 

o.n 

- 1 • 704 ? F-02 

0.0 


2 

0.0 

0.0 

-7 . 868 F- 

-01 

0,0 

-4,nftl6F-02 

0.0 

5 

1 

0.0 

0.0 

— 2 • 1 ^ 3 F 

00 

4 *07 l ^ F -0 4 

-1 . 1643F-02 

0.0 


2 

0.0 

0.0 

-l.?44^ 

00 

l .40 7 8 c — 03 

-4.03?SF-0? 

P.0 

4 

1 

0.0 

0.0 

-2.163F 

00 

0.0 

-l.l670=-02 

0.0 


2 

0.0 

n.o 

-1.244F 

00 

0.0 

-4.0349F-02 

0.0 

3 

l 

n.o 

0.0 

-2.251= 

00 

2.C687F-04 

- 5. 92 7? C -P3 

0.0 


2 

0.0 

0.0 

-1.565S 

no 

7.8716F-04 

-2.2554F-0? 

o.n 

2 

1 

0.0 

0.0 

-2 • 2 = 1 C 

00 

0.0 

- = » 9 30 R F— 0 3 

0.0 


2 

0.0 

o.n 

-1.565F 

no 

o.n 

-2.7567F-02 

0.0 

l 

l 

0.0 

0.0 

-2.281F 

00 

0.0 

0.0 

0.0 


2 

0.0 

0.0 

-L.673F 

00 

0.0 

0.0 

0.0 


***********+**♦**+*** ********** 

EVALUATION OF DESIGN NUMBER 0 

J-****************************’** 


STRESS AREA RATIO 
MAX 0.2392E-03 

MIN 0, ?3fc9F— 03 


LDAP COMP PES VA^TARLF 
1 1 
1 1C 


MAX PUCK RATIO LOAD CCNP 


0.1002E 01 1 

0.H4SF 00 1 


.DESIGN IS CRITICAL 


STRUCTURAL WF I GHT = 0.9692= 02 
REDESIGN OPERATION FOUrws 

OPTIMALITY TNOFX of DESIGN variables FOR RUCKLING CPPS^PAINTS 
DV NO ACT/RAS INDEX 


1 

ACT 

0.17828E 

01 

2 

ACT 

0. 17099F 

01 

3 

ACT 

0. 15695= 

01 

4 

ACT 

0. 13773* 

01 

s 

ACT 

0. 1 1 52 1 F 

01 

ft 

ACT 

0,91532= 

OO 

7 

ACT 

C.ft8834= 

00 

8 

ACT 

0.49002P 

00 

9 

ACT 
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OO 

10 
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00 


P.1.17 
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1 
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ANALYSIS OF DESIGN MJMBF* 5 

******************** A* 4 ******* 


NDQAL DISPLACEMENTS AND ^OTATICNS 


NODE 
NO * 

LOAD 

CASE 

X 

Y 


z 

XX 

YY 

31 

1 

© 

© 

o.c 

0.0 

0.0 

0.0 

0.0 

30 

1 

© 

© 

o.c 

0.0 

0.0 

0.0 

0.0 

29 

1 

© 

© 

0.0 

© 

© 

n.o 

n.o 

0.0 

28 

1 

0.0 

O.o 

0.0 

0.0 

0.0 

0.0 

27 

1 

0.0 

O.c 

0.0 

0.0 

0.0 

0.0 

26 

l 

0.0 

0.0 

© 

o 

n.o 

0.0 

0.0 

25 

1 

o.n 

0.0 

0.0 

o.n 

0.0 

0.0 

24 

1 

o 

* 

© 

0.0 

0.0 

0.0 

0.0 

0.0 

23 

1 

o 

o 

0.0 

© 

o 

0.0 

0.0 

0.0 

2? 

1 

0.0 

o.o 

0.0 

0.0 

o.o 

o 

o 

22 

1 

-1.052F-03 

-3.647E-05 

0.0 

o.n 

0.0 

0.0 

20 

l 

-1.051E-03 

o.c 

0.0 

0.0 

0.0 

o 

o 

19 

1 

-8.227E-04 

-2 • G49E-05 

0.0 

o.n 

o.o 

0.0 

18 

1 

-B. 232E-04 

o.o 

0.0 

o.n 

o.o 

0.0 

17 

1 

-6.820F-04 

-2 • 363E-05 

0.0 

0.0 

0.0 

0.0 

16 

1 

-6.824F-04 

0.0 

0.0 

o.o 

0.0 

0.0 

15 

l 

-5.705F-04 

-1.577F-05 

0.0 

0.0 

0.0 

O.o 

14 

1 

-5.709F-C4 

0.0 

0.0 

0.0 

O.n 

0.0 

13 

1 

-4 .738F-04 

-1.641F-05 

0.0 

o.n 

0.0 

0.0 

12 

1 

—4. 74PE-04 

0.0 

0.0 

o.o 

o.o 

0.0 

11 

l 

-3. B58F- P4 

*-1.3 36F-05 

o 

• 

o 

0.0 

o 

o 

0.0 

10 

l 

-3.860E-04 

o.c 

0.0 

o.n 

0.0 

0.0 

9 

I 

-3.036F-04 

-I.C5CF-05 

0.0 

0.0 

n.o 

0.0 

e 

1 

-3.037E-04 

0.0 

o.o 

0.0 

O.o 

n.o 

7 

1 

- ?• 2 c ?F-04 

-7.772F-06 

o.o 

o.C 

n .0 

n.o 
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1 

-2.253F-04 

0.0 

Q 

O 
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0.0 

5 

I 

-1.494F-04 

-5 * 1 33E-06 

0.0 

o.o 
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0.0 

4 

1 

-1.4S4E-04 

0.0 

0.0 

G.O 

0.0 

o.o 

3 

L 

- 7 • 4 76F-05 

—2 • 536E-06 

0.0 

o.o 

0.0 

0.0 

? 

I 

-7.40OF-O5 

0.0 

0.0 

O.o 

0.0 

0.0 

1 

1 

0.0 

0.0 

o 

o 

0.0 

n.o 

o.o 

VALUES 

OF 

DESIGN VARIABLFS 







l 

2 

3 

4 

5 

6 7 A 

0 

0. 

7222E 00 0. 

7117E 00 0. 

6969F 00 

0 • C 7 3 1 F 00 

0.6406F 00 

0.S982E 0O 0.3436P 00 0.4729F 00 


9 in 

0. 37 7?F 00 0.2371E 00 


ANALYSIS OF PLATE/SHEIL FLUENTS f f.C N S T FA CTCE * 1 


element 

FLFMENT 


LOAD /- 


---f 

^FRANF FHF 

>r f s - 


RFMOP 

^'G/TWTSTTNG 

MpMPNTS 

NUMBFP 

THICKNFS* 

C?ND 

NXX 


NYY 


NXY 

MX X 

myy 

MX Y 

1 

0.7222F 

00 

1 

-0.7666E 

01 

-0.7546P 

01 

0.34 76F-03 

0.0 

0.0 

0.0 

2 

0.7U7E 

00 

1 

-0.7555E 

01 

-0.7474E 

01 

-0.7371F-03 

0.0 

o.n 

0.0 

3 

0.6969E 

00 

1 

-0.75C9E 

01 

-C.7405F 

01 

-0.7579F-05 

0.0 

o.o 

0.0 

4 

0. 673 IF 

00 

l 

-0.7459F 

01 

-0 . 7 2 70 c 

01 

-0.4267F-04 

0.0 

0.0 

0.0 

5 

0.64C6F 

00 

1 

-0.73S6E 

01 

-0.7082F 

01 

-0.4410 F-C 4 

0.0 

0.0 

0.0 

6 

0.5982F 

00 

1 

-C . 73 15 F 

01 

-O.603OF 

01 

-0.4921F-04 

0.0 

o.o 

n.o 

7 

0.5436F 

00 

1 

>0.72146 

01 

-0.6497F 

01 

-0.2819 F— 04 

0.0 

o.o 

o.o 

8 

0.4729E 

00 

1 

-0.7088F 

01 

-0.6043F 

01 

-0 .5957F-04 

0.0 

0.0 

0.0 

9 

0.3772E 

00 

1 

-0.6S25F 

01 

-0.53R7F 

01 

0.1378F-03 

0.0 

O.o 

0.0 

10 

0.2371F 

00 

1 

-G.6706F 

01 

-0.4330E 

01 

-0.4602F-03 

o.o 

O.n 

0.0 


ANALYSIS OF BOUNDARY FLEMEMS - CCNSTRAINT FCBCFS 


CLNST NUMBER 

LOAD C A S p 

Ff?R CF 


1 

l 

0.0 


o.n 

2 

1 

0.0 


0.0 

3 

l 

0.0 


c.o 

4 

1 

0.0 


0.0 

5 

1 

0.0 


0.0 

6 

1 

n.o 


0.0 

7 

l 

n.o 


0.0 

P 

1 

o.n 


c.o 

9 

1 

0.0 


0.0 

10 

1 

0.0 


n.o 

LI 

1 

0 . 750 15 F 

02 

G.O 

12 

1 

0. 744 1 IF 

02 

c.o 

13 

1 

0. 73347E 

02 

n.o 

14 

1 

0.71722F 

02 

o.o 

15 

1 

0.69526F 

02 

c.o 

16 

1 

0. 66 5 C 7F 

0? 

0.0 

17 

1 

0. 626 53 c 

c? 

0.0 


.1 


P.1.20 



18 

1 

0.570R3F 02 

0.0 

19 

1 

0.49457F CZ 

c.c 

20 

1 

0.Z1B71E 02 

0.0 

BUCKLING LOAO 

PAPAMFTFPS 



0.99759O 00 

0.4Z9560 01 



BUCKLING KDOE SHAPES 

NO OF MQDF X 

NO. SHAPE 

31 l 0,0 

2 0,0 
30 I 0.0 

2 0.0 
29 1 0.0 

? 0.0 

2B l 0.0 

? 0.0 
21 1 0.0 

2 0.0 
2.6 I 0.0 

2 0.0 

25 l 0.0 

2 0.0 
24 1 0.0 

2 0.0 

23 I 0.0 

2 0.0 
22 1 0.0 

2 0.0 
21 l 0.0 

2 0.0 
20 1 0.0 

2 0.0 
19 1 0.0 

2 0.0 
18 1 0.0 

2 0.0 
17 1 0.0 

2 0.0 
16 1 0.0 

2 0.0 
15 1 0.0 

2 0.0 
14 l 0.0 

2 0.0 

13 1 0.0 

2 0.0 
12 1 o.n 

2 0.0 
11 1 0.0 

? 0.0 

10 ! 0.0 

2 0.0 
1 0.0 



Y Z 

XX 

YY 


0.0 

0.0 

O.n 

0.0 

0.0 

o.n 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

o.o 

O.o 

0.0 

0.0 

0.0 

o.o 

o.o 

0.0 

0.0 

0.0 

0.0 

c.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

n.o 

0.0 

O.Q 

0.0 

0.0 

o.o 

c.o 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

n.o 

0.0 

0.0 

n.o 

0.0 

o.o 

0.0 

0.0 

c.o 

0.0 

0.0 

0.0 

o.n 

o.n 

0.0 

0.0 

o.n 

o.o 

0.0 

0.0 

0.0 

0.0 

o.o 

n.o 

0.0 

O.Q 

0.0 

0.0 

o.n 

0.0 

0.0 

0.0 

n.o 

o.o 

0.0 

o.n 

0.0 

0.0 

0.0 

0.0 

o.c 

0.0 

o.o 

o.n 

0.0 

0.0 

0.0 

0.0 

o.n 

0.0 

0.0 

0.0 

c.n 

0.0 

0.0 

0.0 

0.0 

-1.7326F-03 

A.9618E-n? 

0.0 

0.0 

o.o 

2.1<>9 7 c -03 

-6.3002F-02 

0.0 

0.0 

0.0 

n.o 

A.Q67nc-nz 

0.0 

o.c 

0.0 

o.o 

-e.^oenF-n? 

0.0 

0.0 

4.574E-C1 

-1 .33605-03 

?.fl?6Z c -02 

c.o 

0.0 

-A • 1 1 OF- 0 1 

1.663OE-04 

-4.7490F-03 

0.0 

o.n 

4.574E-C1 

n.o 

? .R287F-P2 

0.0 

0.0 

-4.11CF-01 

o.n 

-4.7712F-03 

0.0 

o.c 

0.O62E-O1 

-1.07?7F-(H 

? .07 5 1 C ~C2 

o.o 

o.n 

-3.3Z6E-01 

-7.0464F-04 

7.O177C-0? 

0.0 

0.0 

8.062E-G1 

0.0 

3 .076 9F-0 ? 

0.0 

0.0 

-3.326E-C1 

o.n 

2.0 1 90F-P? 

0.0 

0.0 

l.np ?c on 

-8.7790^-04 

7.514 6 F-P? 

0.0 

0.0 

-7.857F-CZ 

-l.C^fl^P-03 

- ? .9 7^5 F-p? 

0.0 

0.0 

1.087F 00 

0.0 

2.5161^-02 

o.n 

o.n 

-7 • 8 56 E- G 2 

0.0 

7.9754 F-Q? 

o.o 

0.0 

1.31 7 E 00 

-7.163? c -04 

2.P51 7E-CZ 

0.0 

0.0 

2.342E-01 

-l . 1183F-n3 

3.202PF-P2 

0.0 

o.o 

l . ■* 1 00 

0.0 

7.05Z9F-0? 

0.0 

o.n 

2.34? C -C1 

0.0 

3 ,?047 c -r2 

0.0 

0.0 

l.^OZF CO 

-5.7516^-04 

1 . 64 75 C- P2 

0.0 

o.O 

5.490 c -01 

-1 ,n^R7 F-03 

*» .o'* 7*F-r? 

0.0 

o.O 

1 • 50? r cn 

O.n 

1 .f 4P5P-0? 

O.o 

0.0 

5.490E-01 

0 .0 

3.0343T-P2 

0.0 

0.0 

1.649P 00 

-4.4 7^nr-n4 


o.n 


9 


n 





2 

o. 0 

o.o 

8 • 342 c - 

01 

-9. 1607F-04 

7.A241F-Q? 

o.n 

B 

1 

o.o 

0.0 

1.649E 

00 

0.0 

1.2P24F-02 

0.0 


2 

o.o 

o.o 

8.342F- 

01 

‘0.0 

2.6257F-C? 

0.0 

7 

1 

0.0 

n.C 

1.761F 

00 

-3.2B7Or-04 

9 , A 1 7 1 f-03 

n.o 


2 

0.0 

o.o 

1.071F 

On 

-7 .2214F-04 

2 .o^pRr-02 

0.0 

ft 

I 

0.0 

o.o 

1. 76 If 

00 

0.0 

9.^2?P f: -03 

0.0 


2 

0.0 

0.0 

1.071*= 

CO 

0.0 

2.0701F-H2 

0.0 

5 

1 

o.o 

o.c 

1.839? 

00 

-2.1M3F-C4 

fc.l933F-03 

O.o 


2 

o.o 

0.0 

l.?*f r - 

00 

-4 .*=A72F-04 

1 .4??*F-0? 

O.o 

9 

1 

0.0 

o.c 

1.R3RF 

00 

o.o 

6.1971F-03 

0.0 


2 

0.0 

0.0 

1.2RFF 

00 

0.0 

1.4242F-02 

o.o 

3 

l 

0.0 

o.o 

1.685*? 

00 

-1 .06ATF-CK 

' a .057?F-n 3 

0.0 


? 

0.0 

O.o 

1.394F 

CO 

-?.snnF-o<v 

7.2O80F-03 

0.0 

2 

1 

0.0 

0.0 

1 .BOSE 

00 

0.0 

3 . 099 2F-02 

0.0 


2 

n.o 

o.n 

1 . 3 5 4 F 

00 

o.n 

7. 212 A f-03 

0.0 

1 

1 

0.0 

0.0 

1.O01E 

no 

0 .0 

C.O 

0.0 


2 

0-0 

o.c 

1 • 39 1 F 

no 

0 .0 

n.o 

0.0 
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EVALUATION of DESIGN NUMBER *> 

*****4******4*<-*4**4-<*-**4*4 4**** j 
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STRESS AREA RATIO 
MAX 0.4969F-03 

MIN 0.2107F-03 

MAX BUCK PATTp 


ICAn C CNC 
1 
1 


CES VAR I API F 
10 

1 


LTAO CTNC 


0 , 1 002F 01 1 

0.2328F O0 1 


CFSIGN IS CRITICAL 


STRUCTURAL WF I GHT= 0.8459F 02 


REDESIGN OPERATION FOLLOWS 


OPTIMALITY INPFX 0= DFSIGN VARIATE FOP PL CK L I AC CONSTRAINTS 
DV NO ACT/PAS INTEX 


1 

ACT 

0.9P7<*4E 

00 

2 

AFT 

C.9BB9BF 

00 

3 

ACT 

C.99151F 

nr 

4 

ACT 

G.99373F 

or 

9 

ACT 

0.9965ZF 

rc 

6 

ACT 

C. 10001F 

01 

7 

ACT 

0. 100 5 IF 

ni 

fl 

AC T 

0. 1 0 1 32 c 

ot 

Q 

ACT 

C. 10295*: 

oi 

10 

ACT 

0 • 9E A 73F 

On 


P.1.22 



NO. OF ACTIVE BUCKLING CONSTRAINTS Apc l 
Bl/CKLING - CRITICAL DESIGN HAS CCAVFRGfr 
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